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ABSTRACT 

All  present  stiffened  cylindrical  shell  design  formulas 
for  the  case  of  external  hydrostatic  pressure  were  surveyed. 
This  also  included  present  design  practices  for  allowances 
due  to  pressure  hull  imperfections,  and  actual  test  data 
when  available.   Formulations  for  all. three  basic  hull 
failure  modes  were  then  selected ,  first  for  accuracy,  and 
secondly  for  compatibility  with  both  elastic-perf ectly 
plastic  and  strain-hardening  metals,  whenever  possible. 

The  formulas  were  then  inserted  in  a  logical  flow 
pattern  to  design  elastic-perf ectly  plastic  scantlings  .for 
failure  at  the  most  efficient  collapse  mode.   The  process 
was  programmed  in  FORTRAN  IV  and  designed  to  iterate,  varying 
several  scantling  parameters  systeio.atica.lly.   The 
''optimum"  design,  based  on  a  simple  hull  weight/buoyancy 


ratio  was  selected  at  the  completion  of  the  run. 

Program  inputs  are:   collapse  depth,  hull  diameter, 
hull  length,  internal  bulkhead  spacing  (specified  or 
unspecified),  framing  (internal  or  external),  and  metal 
properties.   Outputs  are:   shell  thickness,  typical  frame 
spacing,  typical  frame  size,  heavy  frame  (bulkhead)  spacing, 
and  heavy  frame  si^e  for  each  design,  and  an  optimum  design 
designation.   Simple  directions  are  given  for  conversion 
of  the  program  to  one  compatible  with  strain-hardening  metals, 
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INTRODUCTION 

The  amount  of  literature  concerning  the  collapse  of  ring- 
stiffened  cylindrical  shells  under  hydrostatic  pressure 
accumulated  in  the  last  fifty  years  is  voluminous.   This  is 
understandable;  the  subject  is  very  involved.   To  this  day, 
an  exact  solution  for  all  aspects  of  the  problem  does  not 
exist.   Good  solutions  for  the  different  failure  modes  do 
exist,  though,  modified  in  varying  amounts  by  empirical  data. 
No  attempt  will  be  made  here  to  list  or  summarize  this 
knowledge.   liany  have  already  done  this,  and  the  finest 
review  to  date  in  this  author's  knowledge  has  been  done  by 
J.  G.  Pulos   for  the  Navy's  former  David  Taylor  Model  Basin. 

It  appeared  that  one  should  be  able  to  integrate  this 
albeit  incomplete,  yet  extensive  knowledge  with  the  use  of 
present  generation  computer  science.   Hand  calculations  for 
only  one  geometrical  combination  of  shell  thickness,  frame 
size,  frame  spacing,  hull  diameter,  hull  length,  etc.  are 
notoriously  laborious,  even  for  only  one  mode  of  failure. 
Submarine  design  processes  using  the  hand  technique  'would 
achieve  adequate  structures,  but  with  little  or  no  idea  if 
anything  better  existed. ■  Optimization,  with  the  exception  of 
a  few  combinations  tried  at  a  great  cost  in  time,  was  out  of 
the  question:   while  similar  submarines  could  be  designed  on 
past  knowledge,  different  hull  geometries  or  deeper  operating- 
depths  meant  a  great  deal  of  time  and  work.   It  was  at  the 


16 
suggestion  of  Professor  Evans  that  the  development  of  a 
computerized  design  optimization  was  undertaken.   The  basic 
design  equations  were  there.   The  computing  tools  were 
readily  available.   All  that  needed  to  be  done  was  an 
integration  of  the  two  into  a  practical,  useable,  and  most  of 
all i  reliable  (in  terms  of  latest  empirical  data,  if 
necessary)  program. 

Shortly  after  the  start  of  the  project,  it  was 

op 

discovered  that  a  smiliar  program  had  just  been  completed^ 
for  the  Naval  Ships  Research  and  Development  Center  (NSRDC). 
It  is  hoped  that  by  using  some  different  approaches  and 
techniques,  this  program  could  be  a  valuable  tool  to  use  in 
conjunction  with  reference  22. 

The  description  of  the  program  development  will  be  done 
by  subprograms,  each  building  on  the  other,  and  ending  with 
the  optimization  scheme  of  the  main  program.   While  the 
program  developed  can  be  used  only  with  elastic-perf ectly 
plastic,  isotropic  materials  (e.g.,  HY-80  steel),  it  can 
easily  be  modified  so  as  to  be  applicable  also  to  strain- 
hardening  metals,  such  as  HY-150  steel.   Further  discussion 
on  this  will  follow  later,,   Included  also  in  the  thesis  will 
be  various  data  obtained  using  the  program  or  portions 
thereof  in  parametric -type  studies « 
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It  Is  recommended  that  while  reading  through  the  various 
chapters  on  subprograms,  reference  be  made  to  figure  23 
(chapter  10),  which  is  a  simplified  main  flow  diagram  for 
the  entire  optimization. 


18 

CHAPTEE  ONE 
BACKGROUND  INFOP.riA.T10N  AND  TERMINOLOGY 
The  entire  hull  design  program  is  based  on  the  three 
fundamental  failure  modes  for  ring-stiffened  cylindrical 
shells.   Therefore,  in  order  that  terminology  remain 
consistently  clear  throughout  the  discussion,  a  brief 
description  and  catagorization  of  these  modes  follows.   Refer 
to  figure  one  for  pictorial  representations, 

A«   AXI SYMMETRIC,  FAILURES BETWEEN  STIFFEKEES 

All  axlsymmetric  failures,  whether  elastic,  plastic, 
or  some  mixture  thereof,  are  characterized  by  one  or  more 
accordian-like  pleats,  or  circumferential  ripples  between 
ring  frames.   For  true  axi symmetric  failures,  the  stiff eners 
remain  undeformed. 

!•  AXISYKKSTRIC  YIELD  FAILURE.   This  type  of  failure 
occurs  only  with  elastic-perf ectly  plastic 
(plateau-type  stress-strain  curve)  materials „   It 
is  regarded,  then,  as  almost  totally  a  yield-type 
failure,  although  it  is  initiated  partially  by 
instability  phenomena. 
2«  AXI  SYMMETRIC  -I]     STIC  FAILURE.   Inelastic 
failure  is  also  a  failure  above  the  purely 
elastic  range,  but  in  strain-hardening  materials. 
Thus,  the  failure  is  in  the  range  where  Young's 
modulus  varies ,  and  can  intuitively  be 


considered  a  kind  of  combination  elastic-plastic 

failure.   Often  this  failure  can  be  at  a  lower 

pressure  than  a  pure  buckling  (elastic)  failure, 

d.ue  to  certain  combinations  of  lov?  modulus  and 

hull  geometry.   The  pure  yield  failure  (A.l.) 

cannot  occur.  Ah  strain-hardening  materials. 

3*  AXI SYMMETRIC  ELASTIC  FAILURE.   This  is 

ax i symmetric  failure  in  the  linearly  elastic 

range.   Theoretically,  this  could  occur,  given 

proper  geome.try,  in  either  of  the  two  above  types 

of  materials.   Generally  speaking,  the  required 

geometry  is  one  of  a  thin  shell  relative  to  the 

hull  diameter  and  depth.   In  reality,  this 

failure  is,  at  this  writing,  a  mathematical 

phenomenon  only.   Other  hull  failure  modes  occur 

first,  so  this  mode  has  never  been  achieved  in 

actual  testing.   It  is  valuable,  however,  in 

determining  effects  of  geometrical  defects  on 

collapse  pressure. 

B .  ASYM M Ejny^^LCBAR )  BUCKLING  FAILURES  BETWEEN 
STIF^NERS: 

All  lobar  buckling  failures  ore  characterized  by- 
lobes  of  buckling  distributed  partly  or  completely  around  the 
shell  circumference.   The  frames  remain  intact. 

1.  ELASTIC  OR-  INE31ASTIC  HOPES:   This  is  primarily  a 
buckling  failure.   The  collapse  pressure  is 
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dictated  by  the  modulus  (elastic  or  inelastic) 
and  hull  geometry, 

c.  gf^ier/l  instability _,(sk5ll  a^d_f?:aks)..,faii,ufe. 

General  instability  failures  are  characterized  by 
failure  of  both  frames  and  shell  simultaneously,  sometimes 
extending  the  entire  length  of  the  cylinder. 

1.  ELASTIC  OP.  INELASTIC  HOPES.   Ihese  are  also 
,   buckling  failures,  except  that  lobes  extend 
longitudinally  as  well  as  circumf erontially. 
It  is  normally  assumed  that  only  a  one-half  wave 
extends  between  heavy  frames  (bulkheads),  since 
these  heavy  frames  are  designed  heavy  enough  not 
to  deform  at  the  general  instability  collapse 
pressure.   This  combination  failure  of  frames 
and  shell  is  not  as  well  understood  or  formulated 
mathematically  as  the  other  two  modes,  especially 
in  the  inelastic  region.   Thus,  larger  safety 
factors  are  employed  when  checking  hull  designs 
in  this  mode. 
It  should  be  emphasized  at  the  outset  that  the  program 
to  be  described  designs  t?ie  stiffened  cylindrical  shells  to 
fail  in  the  ax i symmetric  yield  failure  mode.   For  the 
untested  inelastic  portion  of  the  program,  the  design  failure 
mode  would  be  axisymmetric  inelastic  failure B   This 
philosophy  is  used  today  in  submarine  design  and  is 
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-advocated  by  most  naval  archietects.   It  has  tv?o  basic 
reasons: 

1)  Failure  by  yielding  of  the  shell  utilizes  full 
material  (i.e.,  yield)  strength.   Stresses  in 
buckling  type  failures  are  usually  below  yield 
stress,  -with  failure  depending  mainly  on  the  value 
of  Young's  modulus  and  geometry  of  the  failing 
structure.   More  efficient  structures  should  thus 
result  from  designing  to  a  yielding  failure. 

2)  Imperfections  in  construction  (e.go,  hull  out-of- 

roundness)  effect  buckling  failures  far  more 

? 

seriously  than  yield  failures.    Thus,  hulls 

designed  for  yield  failure  would  have  less 
stringent  requirements  for  building,  and  less 
chance  of  failure,  given  that  imperfections 
might  exist. 

The  optimization  criterion  used  is  a  simple 
hull  weight  to  hull  displacement  ratio,  which 
appears  to  be  the  best  general  measure  of  design 
efficiency  for  this  type  of  structure. 


Figure  1. 
FAILURE  MODES:  STIFFENED  CYLINDRICAL  SHELLS 
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AXISYMMETRIC  YIELDING  (BETWEEN  FRAMES) 
COLLAPSE  FORMULAS:  LUNCHIK ( DTMB#  1291.1393) 


ASYNHSTRI C ( LOBAR ) BUCKLING ( BETWEEN  FRAMES ) 
COLLAPSE  FORMULA:  REYNOLDS (DTMB#1 392 ) 


GENERAL  INSTABILITY  (CONCURRENT  FAILURE 
OF  SHELL  AND  FRAMES) 

COLLAPSE  FORMULA:  KHENZKE-KIERNAN-(DTMB#  16??) 
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CHAPTER  TWO 
HULL  STRESSES 
Numerous  subprograms  throughout  the  optimization  require 
the  values  of  various  stresses  in  the  hull.   These  stresses 
include  circumferential  and  longitudinal  stresses,  at  the 
frames  and  at  midbay ,  at  the  inner  and  outer  shell  surfaces. 
Several  solutions  for  stresses  due  to  external  hydrostatic 
pressure  on  ring-stiffened  cylinders  have  appeared  in  the 
past.   The  most  famous  was  due  to  the  Germans,  Von  Sanden 
and  Gunther,  in  1920.   Portions  of  their  analysis  are  still 
in  use  today.   In  1930 •  the  Italian  Viterbo  modified  their 
analysis  to  include  the  so-called  stiff ener-expansion  effect 
(a  result  of  axial  stresses  in  the  shell).   Neither  of 
these  early  analyses,  however,  included  the  "beam  column" 
effect.   This  effect,  introduced  by  Salerno  and  Pules   in 
1951,  is  caused  by  the  interaction  of  longitudinal  bending 
and  longitudinal  compression  in  the  hull  caused  by  the  axial 
portion  of  the  hydrostatic  pressure  acting  on  it.   The 
Salerno-Pulos  stresses  are  an  exact  solution,  and  the  beam- 
column  effect  accounts  for  any  non-linearities  between, 
pressure  and  strain  in  the  cylinders.   In  all  cases  except 
one  (see  chapter  four),  the  program  uses  the  more  accurate 
Salerno-Pulos  (hereafter  S-P)  stresses. 


26 
The  beam- column  effect  is  represented,  in  the  S-P 
analysis  by  the  parameters ( hereafter:   GAMA),  where  : 


^  =  -!r*=£i/3(7^  (t1) 


p*  is  defined"  '"  '  as  the  "critical  load  for 
axi symmetric  elastic  buckling  of  an  unstiffened  cylindrical 
shell  under  the  action  of  uniform  axial  pressure. "  GAMA=0 
corresponds  to  a  zero  beam-column  effect  (i.e.,  the  stress 
solution  of  Von  Sanden  and  Gunther,  hereafter  V-G).   As 
GAMA  grows  larger,  the  beam -column  effect,  and  thus  the  non- 
linearity  betv.'een  pressure  and  hull  stresses,  increases. 
When  GAMA^>  1,  then  theoretically  the  between-frame  failure 
mode  shifts  to  axi symmetric  elastic  buckling  (see  chapter 
one  for  definition).   As  explained  in  Chapter  1,  this  appears 
to  be  a  mathematical  phenomenon  only,  for  it  has  never  been 
achieved  in  actual  testing*   Other  modes  of  failure  (e.g., 

lobar  buckling)  always  appear  to  occur  first,  or  else  the 

'i  3 

axisymmetric  failure  is  always  accompained  by  some  yielding. 

The  S-?  stresses  are  calculated  by  placing  various 
combinations  of  S-P  "F  functions"  (nee  reference  18  for 
expressions  and  curves)  into  the  S-P  stress  equations. 
Formerly,  this  was  a  very  laborious  process,  far  more  time- 
consuming  than  obtaining  stresses  through  the  simpler  V-G 
equations.   The  curves  developed  by  M.A.  Krenzke  and  E.  D. 
Short  (included  in  reference  18)  shortened  the  labor 


2? 


considerably.   The  computerized  solution  makes  it  almost 
mandatory  to  use  the  superior  S-P  stress  solution. 

The  "F  functions"  are  obtained  through  the  use  of 
two  basic  subprograms  (see  figures  2,    3  and  h) t      Subroutine 
FRAME  essentially  takes  shell  thickness,  unsupported  shell 
length,  and  frame  moment  of  inertia  and  manipulates  them  to 
obtain  the  S-P  variables  THETA,  ALFA,  and  BETA  (see  Appendix 
A  for  method  of  obtaining  frame  dimensions).   These  variables 
are  then  transferred  to  .the  stress  program  PULOS  via  a 
COMMON  statement.   FRAME  is  separate  from  PULOS  because  in 
one  program,  (HVYFRM)  PULOS  is  used  in  two  different  places 

with  the  same  scantlings.   All  variables  in  FRAME  are 

18 

computed  in  accordance  with  the  S-P  stress  analysis  v ,  with 

the  exception  of  the  expression  for  effective  area,  AE. 
Reference  18  lists  this  as  (in  program  terms): 

AE=AF*(RM/RF),  or 

AE-AF •::•  ( E  M/R  F )  *  ■'-:■  2  ,  d  e  p e nd  1  ng  o n  wh ether 

the  framing  is  internal  or  external,  respectively.   This  is 

oh 

not  strictly  correct,  and  has  since  been  refined  by  Short'  , 

who  used  the  similar  equation: 

AE=AF#(RM/RF)**Q,  where  Q=l+  2*GKU 
This  equation  is  good  for  either  internally  or  externally 
framed  cylinders  with  "reasonable"  (i.e.  ,"  suitable  for  this 
program's  purposes)  frame  depth/shell  radius  ratios. 

The  subroutine  PULOS  is  a  straightforward  adapation 
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of  the  S-P  stress  analysis.   It  produces  the  four  "F 

functions11  and  the  variable  "a"  (see  reference  16),  the 

various  combinations  of  which  are  used  in  several  parent 

programs  to . compute  the  hull  stresses. 

From  the  expression  for  GAMA ,  it  can  be  seen  that  an 
input  pressure  is  required  for  PULOS  to  compute  its  outputs. 
In  other  words,  FRAME  and  PULOS  can  be  used  to  calculate 
the  hull  stresses,  given  scantlings,  material  properties, 
and  hydrostatic  pressure.   More  often,  however,  the  program 
is  attempting  to  find  a  critical  failure  pressure.   This 
results  In  two  unknowns,  the  pressure  and  the  stresses,, 
GAMA  thus  becomes  the  result,  of  a  transcendental  process, 
in  which  PULCS  calculates  stresses,  and  the  pressure  input 
is  from  a  parent  program  dealing  with  a  particular  hull 
failure  mode. 

Great  difficulties  were  experienced  in  this  iterative 
process  when  GAMA^l.O;  this  meant  that  ETA1  (see  figures 
3  or  M  became  imaginary  (i.e.,  the  failure  mode  had  shifted 
to  axisymmetric  elastic).   To  combat  this,  several  methods 
were  employed,  each  one  being  different  for  different 
parent  programs.   This  is  the  reason  for  PULCS  and  PULOS  1. 
The  methods  of  validly  circumventing  this  pitfall  will  be 
explained  separately  under  each  individual  failure  mode 
subprogram.  v     ■ 
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CHAPTER  THREE 
COLLAPSE  PRESSURE  REDUCTION 

Realistically  speaking,  no  cylinder  that  is  manufactured 
today  can  be  considered  "perfect".   There  will  always  be 
some  reduction  in  the  collapse  pressure  due  to  manufacturing 
imperfections,  such  as  shell  or  frame  out-of -roundness,  and 
to  residual  stresses  from  welding.   If,  however,  a  pressure 
hull  is  machined  rather  than  rolled  and  welded,  a  structure 
that  is  "perfect"  for  all  practical  purposes  may  be  attained. 
Of  course,  for  large  pressure  hulls,  the  expense  (or  even 
impossibility)  incurred  due  to  size  prohibits  construction 
by  machining  only.   Thus,  some  allowances  must  be  made  in 
scantling  computation. 

The  best  overall  method  (i.e.,  including  perfectly 
plastic  and  strain  hardening  materials)  yet  devised  is 
presented  as  a  graph  in  reference  10  (see  figure  5)»   Here, 
the  lower  curve  represents  an  envelope  of  numerous  model 
test  results  conducted  over  the  years  at  the  Model  Basin. 
It  may  be  observed  that  "the  factor  which  is  all  important 
in  determining  imperfection  sensitivity  is  the  margin  of 
stability  Pe/Pj[i  the  absci!  -  t  on  the  plot.   The  lower  the 
margin  of  stability,  the  greater  the  sensitivity  to 
imperfections".21-   The  reduction  factor  (FEDFAC,  or  pEXP/pExP 
for  machined  models  )  located  on  the  ordinate  is  applied  to 


Figure  5.  -^ 

A  COMPARISON  OF  MODEL  BASIN  DATA  ON  MACHINED  RING- 
STIFFENED  CYLINDERS  WITH  WELDED  CYLINDERS  .  .  . 
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collapse  pressures  of  all  modes  in  the  optimization 

program  (see  figure  6). 

The  envelope  curve  is  approximated  by  the  equation: 
REDFAC=0.6**667+Q.11367*BTO-0.00967*RTO**2 ,  where 

RTO=Pe/Pi 

This  equation  and  the  general  form  of  F.EDPR  were  adopted 

from  a  similar  program  in  reference  22. 

As  mentioned  above,  this  safety  reduction  factor  is 

applicable  not  only  to  this  optimization,  but  also  to  one 

for  strain  hardening  metals  as  well.   For  this  reason  (i.e., 

the  fact  it  is  easily  convertible  to  the  strain  hardening 

case)  it  was  chosen  above  other  existing  out-of -roundness 

analyses  for  the  strictly  elastic  buckling  cases,   No 

mathematical  analysis  presently  exists  for  the  inelastic 

21 

case,  due  to  its  complexity.    Further  discussion  on  this 
is  contained  in  Chapter  13. 
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CHAPTER  FOUR 
ASYMMETRIC  (LOBAR)  BUCKLING 

The  most  widely  used  formulation  for  this  failure  mode 
over  the  years  has  been  the  so-called  "DTMB  Instability 
Formula",  developed  by  Vindenburg  and  Trilling.   The  formula 
was  good  only  for  an  elastically-perf ectly  plastic  material. 
More  recently,  Reynolds  v  has  developed  a  more  generalized 
formulation  which  may  be  used  for  either  elastic-perfectly 
plastic  or  strain-hardening  materials.   In  reference  20, 
Reynolds  recommended  using  the  V~G  stresses  (longitudinal 
and  circumferential  at  mid-bay).   He  stated  that  the  accuracy 
of  the  analysis  was  not  seriously  impaired  by  not  using  the 
more  cumbersome,  yet  more  accurate  3-P  stresses.   This  is 
graphically  borne  out  by  the  example  given  in  figure  7» 
In  any  case,  the  theory  correlates  very  closely  with 
experiment,  to  within  four  percent  . 

The  S~P  stress  programs  (FRAME  and  PULOS)  can  easily  be 
made  common  with  any  other  subprogram.   Originally  it  was 
decided  that  there  would  be  redundancy  involved  if  the  less 
accurate  V-G  stresses  were  used  for  the  Reynolds  lobar 
buckling  subprogram  (RNLDS).   Difficulties,  however,  were 
immediately  encountered  when  the  S-P  stresses  were  utilized. 
RNLDS  is  used  in  the  main  program  chiefly  as  a  checking 
function..   Since  the  hull  scantlings,  as  mentioned  in 
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.chapter  one,  are  designed  for  yield  failure,  RNLDS  is  used 
only  to  insure  that  the  scantling  set  under  scrutiny  does 
not  fail  by  lobar  buckling.   Thus,  since  scantling  sets 
being  examined  will  usually  fail  first  by  axi symmetric 
yielding,  failure  by  lobar  buckling  sometimes  occurs  at  much 
greater  pressures.   In  many  cases,  failure  by  axisymmetric 
e la st i c  buc  kl i ng  will  occur  at  a  pressure  between  failure  by 
yield  and  failure  by  lobar  buckling.   As  noted  in  chapter 
two,  this  phenomenon  causes  GAI-1A  (beam-column  effect)  to 
increase  in  value  over  1.0,  resulting  in  imaginary  values 
occuring  vjithin  the  PULOS  subprograms.   Efforts  to  "force" 
convergence  of  RNLDS  by  inserting,  for  instance,  values  of 
GA11A=0.95.  or  of  taking  only  absolute  values  of  the  radical 

/l.0«GAi!A,  were  only  mildly  successful.   The  value  of  final 
convergence  in  any  ca.se  was  not  accurate,  and  certainly  was 
not  that  of  lobar  collapse  pressure. 

For  these  reasons,  then,  V-G  stresses  were  used  in 
RNLDS  (see  figure  8),  as  originally  recommended  in  reference 
20.   One  distinct  advantage  of  the  V-G  stresses  is  that 
they  require  no  iteration  for  convergence.   There  is  no 
separate  stress  program  needed ,  and  the  V-G  stresses  are 
directly  (and  quickly)  computed  within  RNLDS.   Since  (see 
chapter  ten)  RNLDS  is  used  itself  in  an  iterative  process 
within  the  main  orograto,  this  results  in  substantial  savings 
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in  computer  time. 

Function  RNLDS,  in  this  case  for  elastic-perfectly 
plastic  material,  predicts  strictly  an  elastic  type  failure. 
Thus,  the  full  pressure  reduction  from  REDPR  (see  chapter 
three)  is  employed  directly  to  compensate  for  imperfections 
and  residual  stresses. 
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.CHAPTER  FIVE 
AXI SYIU  lETRl  C  _(  Y I  ELD )  ^FAILURE 

The  real  core  of  the  optimization  program  is  the 
axisymmetric  failure  mode,  for  the  program  designs  its 
scantlings  to  fail  by  yield  (see  chapter  one),   Three 
subprograms  in  addition  to  FRAME,  PULOS,  and  REDPR  are 
included  in  this  grouping:   PLNCK,  ELNCK  and  THKNS  (see 
figures  9,  10,  and  11). 

The  rather  famous  Von  Sander]  and  Gunther  formulas  92 
and  92A  (utilizing  the  maximum  shell  stress  theory  of 
Rankine  at  the  frame  and  midbay,  respectively)  were  used  in 
design  for  many  years.   Recently,  however,  many  more  solutions 
have  appeared.   With  the  advent  of  S-P  stresses1",  the 
stress  analysis  alone  has  improved  in  accuracy.   The  manner 
in  which  the  stresses  are  used  to  predict  collapse  by 
axisymmetric  yield  varies  greatly.   Generally  speaking,  the 
maximum  strain  energy  theory  of  Mises  and  Hencky  provides 
the  best  manner  of  stress  combination  to  predict  failures 
within  the  various  shell  yield  formulations.   The  point  at 
which  this  is  applied  is  also  subject  to  discussion. 
Although  it  is  generally  agreed  that  the  largest  stresses 
actually  occur  at  .the  frames,  it  is  becoming  evident  that 
data  indicate  the  best  predictors  use  the  mid-bay  area  as 
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the  initiation  point  of  yielding  failures.   The  most 
frequently  accurate  analysis  for  axisymmetric  yielding 
failure  is  that  due  to  Lunchik  "  ' .   Although  his  formulation 
has  not  been  tested  through  complete  ranges  of  geometries 
and  depths,  the  tests  that  have  been  made  indicate  his 
solutions  are  at  least  as  accurate  as  any  others  (to  within 
1%   of  actual  failure  pressure  in  many  cases),  and  much  better 
than  92  or  92A.   In  reference  1^!-,  Lunchik  shows  that  very 
successful  correlations  were  obtained  with  tests  of  ring- 
stiffened  cylinders  ranging  from  A  (thinness  ratio) =0,^-1  to 
A =0.70.   He  recommended  his  formulation,  however,  only  for 
•'cylinders  where  geometries  are  in  the  range  of  axisymmetric 
yielding",  precisely  the  case  in  this  program.   Basically, 
Lunchik  assumes  a  standard  three-hinge  failure  mechanism, 
postulating  that  the  frame  plastic  hieges  fail  first,  and 
predicting  the  pressure  at  which  the  mid-bay  hinge  is 
complete.   The  basic  difference  between  Lunchik' s  analysis 
and  others  is  his  computation  of  "plastic  reserve  strength". 
Plis  structure  does  not  fail  when  some  outer  hull  fiber  at 
raid-bay  has  reached  yield  stress.   It  fails  only  after. this 
plasticity  has  progressed  through  the  shell  at  that  point 
far  enough  to  produce  a  hinge  and  precipitate  failure. 

FUNCTION  VL  ;G.  The  use  of  Function  REPPfl. (see  chapter 
three)  requires  both  an  elastic  and  an  inelastic  collapse 
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pressure  for  the  particular  failure  mode  under  examination. 
For  the  elastic-perfectly  plastic  materials,  there  is  no 
"inelastic  failure  pressure"  as  such,  since  no   strain- 
hardening  is  involved.   Thus,  the  yield  failure  pressure 
outlined  above  is  substituted.   For  the  solution  of  the 
elasti_c  axi symmetric  failure  pressure,  there  remain  two 
possibilities.   One  is  the  "exact"  solution  offered  by  the 
S-P  stress  analysis.   When  GAMA^l.O  the  elastic  ax  i'sy  nine  trie 
mode  occurs.   However,  the  solution  for  this  pressure  is 
bound  up  in  a  new  "F  function"  and  requires  a  rather 
involved  iteration.   Because  the  REDPE  method  is  approximate 
in  any  case,  e.n   exact  elastic  axisymmetric  value  is  not 
required.   Thus,  the  S~?  solution  was  rejected  in  favor  of 
Lunchil'^s  inelastic  axisymmetric  analysis  -'.   This  analysis 
is  very  similar  to  Reynolds'  analysis,  for  asymmetric 
buckling  (see  chapter  four),  in  that  it  can  be  applied  in 
the  strain  hardening  (inelastic  case)  using  the  secant  and 
tangent  moduli.   By  setting  ETAH=0  and  ESEC=E,  the  solution 
breaks  down  to  one  for  elastic  axisymmetric  buckling.   This 
is  the  pressure  computed  by  Function  PLNCK,  and  the  process 
is  taken  directly  from  reference  13* 

SUB? OUTI BE  EL; i Ck   The  "inelastic"  pressure  used  in 
REDPE  (and  also  the  pressure  to  be  reduced  itself)  is 
computed  in  this  subroutine.   As  explained  above,  Lunchik's 
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yield  analysis  for  elastic-perf ectly  plastic  materials 
is  used.   S-P  stresses  are  used  in  the  formulation  (see 
figure  10).   This,  of  course,  meant  an  iterative  process, 
and  occasionally  convergence  problems  were  encountered. 
If  input  geometries  were  satisfactory,  convergence  was 
accomplished  in  four  or  five  cycles.   Occasionally,  the 
iterative  program  (THKNS)  calling  ELNCK  jumped  outside  the 
range  of  convergent  geometries  in  its  search  process.   For 
such  cases,  an  iteration  limit  of  10  was  put  into  ELNCK. 
This  meant  that  the  pressure  going  back  to  THKNS  was  not 
entirely  accurate,  but  good  enough  to  continue  with  a  search 
pattern  to  find  a  convergent  geometry.   An  additional 
"safety  valve"  was  built  into  Subroutine  PULOS  (see  chapter 
two)  to  prevent  GAHA^l.O  and  thus  producing  imaginary  values. 
If  GAHA>1.0,  GAIIA  was  set  equal  to  0.95,  a"d  the  shell 
thjLcJmejis  adjusted  to  achieve  this.   Thus,  depending  on 
entering  geometries,  occasionally  shell  thickness  itself  is 
adjusted  within  ELNCK.   This  is  justifiable,  in  that  the 
overall  program  is  designing  to  a  yield,  not  a  buckling 
failure.   Any  buckling  geometry,  even  if  it  could  be 
converged  upon,  would  not  be  desired.   Lunchik's  "yield 
pressure",  PY  (reference  16,  equation  35)  is  the  result  of 
the  S~P  stresses  obtained,  and  is  the  pressure  at  which 
yielding  begins  at  the  outer  hull  fiber,  at  mid-bay.   His 
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"plastic  reserve  strength  ratio,  PCLE1/PY  (FCTB)  is  then 
computed  and  multiplied  by  PY  to  obtain  Lunchik's  final 
collapse  pressure,  PCLE1.   The  arrival  at  certain  of 
Lunchik's  equations,  made  somewhat  confusing  by  a  misprint 
in  reference  lh- ,    is  done  in  detail  in  Appendix  B.   Before 
sending  the  collapse  pressure  to  THKNS,  it  is  reduced  for 
residual  stresses  and  manufacturing  defects  by  REDPR, 

At  this  point  (see  figure  12),  it  is  interesting  to 
see,  at  least  in  one  case,  how  the  Lunchik  and  Reynolds 
analyses  compare  with  the  log-log  plot  of  hoop  stress  vs. 
the  Windenburg-- Trilling  formula  presented  in  reference  8. 
For  the  particular  hull  diameter  chosen,  Reynolds1  pressures 
follow  Windenburg's  almost  exactly.   Lunchik's  pressures 
show  hoop  stress,  at  least  in  this  case,  to  be  rather 
conservative. 

SUBROUTINE  . THKNS   This  subroutine  uses  ELNCK  in  an 
.iterative  process  to  converge  on  an  exact  shell  thickness 
which  will  fail  by  axisyrametric  yield  at  the  desired 
collapse  pressure.   If  THKNS  cannot  converge  on  a  thickness, 
it  is  obvious  that  the  input  scantlings,  despite  changes 
in  thickness,  are  such  that  failure  by  elastic  axisymmetric  • 
■buckling  occurs  before  axisymmetric  yield  failure.   Since 
the  full  strength  of  this  metal  is  not  then  being  used, 
this  type  of  solution  obviously  is  not  desired.   In  such 
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a  case,  THKNS  transfers  control  to  the  next  iteration  loop 
of  the  main  program  (see  chapter  ten). 
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CHAPTER  SIX 
HULL  THICKNESS  INCREASE 

Generally  speaking,  the  safety  factor  introduced  by 
REDPR  (see  chapter  three)  to 'compensate  for  "imperfections 
and  residual  stresses"  will  be  adequate.   However,  it  should 
be  remembered  that  the  collapse  pressures  developed  thus 
far  all  'are  "triggered"  by  stress  values  at  raidbay.   The 
hJLghest  stresses  actually  encountered  usually  occur  at  the 
frame  faying  "flange".   To  account  for  this,  stresses  in  that 
area  are  limited  to  75.j   of  yield  stress  at  operating  depth 
(assumed  here  to  be  2/3  collapse  depth)  by  subroutine 
STRTHK,  ss  recommended  in  reference  10.   This  is  done  to 
account  for  such  things  as  lovi-cycle  fatigue,  creep,  stress 
corrosion,  and  to  insure  a  reasonable  stress  level  in  the 
frame  flange  prior  to  collapse  for  those  frames  with  an 
initial  out-of -roundness   . 

Subroutine  STRTHK  (see  figure  13)  uses  the  S-P  stress 

in 
analysis    to  compute  all  four  stresses  of  interest  in  the 

shell  at  the  frames:   inner  and  outer  plate  surfaces,  and 

longitudinal  and  circumferential  stress  directions.   If  the 

largest  of  these  is  greater  than  7 5A   yield  stress  at  2/3 

collapse  pressure,  shell  thickness  is  increased  in 

increments  of  $%   until  the  criterion  is  met. 
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CHAPTER  SEVEN 
GENERAL  IHSTA B ILITY 

The  problem  of  finding  accurate  collapse  pressures  in 
the  general  instability  node  has  generally  given  analysts 
more  trouble  than  the  first  two  failure  nodes.   The  most 
accurate  elastic  general  instability  failure  analysis  was 
done  by  S,  Kendrick  in  1953  at  Britain's  Naval  Construction 

Research  Establishment,  and  is  generally  known  as  the 

3 

"Kendrick  Part  III"  solution'.   One  year  later,  A.  P.  Bryant, 

working  in  the  same  establishment,  developed  a  far  simpler 

n 

approach'  (Kendrick's  was  exceedingly  complex).   Bryant's 
solution  was  a  single  two-term  equation  which  could  be  (and 
has  been)  used  for  design  studies  without  extensive 
computerization,  although  its  accuracy  left  something  to  be 
desired  (it  was  non-conservative).   Basically  speaking, 
Bryant's  equation  incorporates  the  "split  rigidities 
concept",  where  one  term  represents  the  contribution  of  the 
shell,  and  the  other  the  contribution  of  the  frames  and  a 
frame-length  of  shell  plating.   Although  Kendrick' s  solution 
was  put  in  a  simple  graphical  form  by  Reynolds  J    (later 
extended  by  Ball^) ,  it  was  rejected  for  use  in  this  program 
for  two  reasons: 

1)  since  the  general  instability  pressure  is  used 
merely  to  check  the  solution  designed  for  yield 


5^ 

failure,  extreme  accuracy  was  not  required;  and 
2)  Kendrick's  solution  is  good  only  for  elastic 
failure,  and  cannot  be  applied  to  strain- 
hardening  materials.   Since  this  program  is 
designed  to  be  easily  converted  to  one  x-?hich 
can  also  handle  strain -hardening  materials, 
Kendrick  Part  III  was  rejected. 
Bryant's  solution  was  also  rejected  for  the  second  reason 

above.   More  recently,  a  very  convenient  and  more  accurate 

11 

formula  has  been  developed  by  Krenzke  and  Kiernan  •   It  is 

very  similar  in  structure  to  Bryant's  formula,  but  can  be 
used  for  either  Ideally  plastic  or  strain-hardening  materials. 
At  about  the  same  time,  a  very  similar  solution  was  worked 
out  by  Lunchik  at  the  Model  Basin.   Both  the  Krenzke -Kiernan 
and  the  Lunchik  analyses  give  about  the  sane  results  when 
compared  with  actual  test  data  .   Krenzke -Kiernan  "was 
somewhat  arbitrarily  selected  for  this  program,  solely 
because  it  appeared  to  be  referenced  more  often  in  the 
literature. 

SUBROUTINE  KRIiZK  This  subroutine  is  called  by  the 
general  Instability  subprogram,  and  essentially  computes 
failure  pressure  using  the  Krenzke-Kiernan  formula.   It 
also  returns  the  value  of  "N"  (number  of  circumferential 
collapse  lobes)  which  gives  the  lowest  (most  critical) 
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failure  pressure.   It  might  be  noted  that  this  program,  like 
some  others  in  the  optimization,  is  designed  for  dual  usage: 
by  Ideally  plastic  or  by  strain-hardening  materials  (see 
figure  1^). 

SUBROUTINE  GIKST  (sec  figure  15)  G'INST  tests  the 
incoming  scantling  set  from  the  main  program  for  failure  by 
general  instability.   Depending  on  whether  the  collapse 
pressure1 is  too  shallow  (test  failed)  or  too  deep  (design 
too  conservative),  heavy  frame  (i.e.,  bulkhead)  spacing  is 
adjusted  so  as  to  give  a  collapse  pressure  that  is  either: 

1)  greater  than  1.05  times  desired  collapse  depth 
(this  requirement  for  a  small  safety  margin  is 
due  to  the  uncertainties  of  general  instability 
design,  and  comes  from  reference  10),  or 

2)  less  than  2  times  desired  collapse  depth. 

If,  in  this  process,  bulkhead  spacing  becomes  longer  than 
the  entire  pressure  hull  itself,  BS  is  set  equal  to  the 
hull  length  (HULNTH) ,  and  the  requirement  for  heavy  frames 
is  dropped  (see  chapter  eight  on  Subroutine  HVXFRM). 
Incoming  bulkhead  spacing  may  either  be  declared  or 
undeclared,  in  the  input  data.   If  undeclared,  starting 
bulkhead  spacing  is  taken  as  approximately  twice  the  hull 
diameter,  but  always  a  multiple  of  the  small  frame  spacing. 
The  frame  term  in  the  Krenzke-Kiernan  formula  contains, 
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as  one  of  its  variables,  the  "effective  moment  of  inertia", 
EI.   This  is  defined  as  the  combined  moment  of  inertia  of  a 
small  frame  and  its  effective  length,  EL,  of    shell  plating. 
This  differs  slightly  from  the  Bryant  formula,  which  merely 
uses  a  frame _  space  (FS)  of  shell  plating.   The  EL  is 
determined  via  the  S-P  stress  analysis1-  ,  using  the  equation: 

EL=SL*F1+B 
The' use  of  either  EL  or  IS  in  the  formula  makes  little 
difference  in  mos_t  cases;  however,  since  all  the  tools  viere 
handy,  the  EL  was  used  when  Possible.   Problems  of 
convergence,  however,  again  developed  in  the  PULOS  subprogram. 
When  general  instability  pressure  became  very  large  (and  it 
does  for  many  geometries),  the  value  of  GAMA  (see  chaoter 
two)  exceeded  1.0.   Thus,  a  separate  subprogram  (PULOSl) 
was  added,  which  set  EL=FS  when  GAMA^l.O.   This  is  totally 
acceptable,  since  this  case  occurs  only  when  the  general 
instability  collapse  pressure  is  far  too  deep  to  worry  about 
(its  accuracy  is  only  slightly  degraded  anyway).   When 
successive  iterations  of  GINST  expand.  BS  and  bring  this 
collapse  pressure- to  shallower  depths,  EL  can  again  be 
accurately  computed  by  PULOSl. 
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CHAPTE1  EIGHT 
HJ&YXJ2AKES 

Subroutine  HVYFBM  (see  figure  1'6 )  computes  the  minimum 
size  of  heavy  frames  needed  to  insure  the  general  instability 
collapse  pressure  calculated  by  its  preceeding  subprogram 
GINST.   Until  recently,  heavy  frames  were  designed  using  the 
standard  Levy    formula.   It  was  found  after  considerable 

testing,  however,  that  this  formula  often  gave  unsafe 

17 
estimates   .   A  new  formula  was  derived  at  the  Model  Basin 

by  Blumenberg  in  19&5 •  which  agrees  much  more  closely  with 

testSc   Results  from  testing  indicated  that  the  effectiveness 

of  a  particular  size  of  heavy  frames  decreases  as  the 

cylinder  is  lengthened  and  also  that  the  minimum  size  of 

heavy  frames  needed  to  localize  failure  between  heavy  frames 

is  possibly  not  dependent  upon  their  spacing-5 .   Although 

adequate  testing  has  not  yet  been  published  to  positively 

confirm  Blumenberg 's  formula,  initial  results  show  it  is 

better  than  what  was  formerly  used,  and  thus  it  was  put  into 

HVYFRM: 

FIH=  PCG»]   '  C**3 where- 

tltt      ((M»«2-1.0)*ETAN)'  wnere' 

1)  E  may  be  substituted  for  ETAN  with  ideally  plastic 

materials.   ETAN  actually  is  not  in  Blu  ienbergfs 

formula,  but  was  placed  there  by  the  author  to 

make  HVYFEM  useful  in  an  inelastic  optimization, 
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•  2)  M  is  the  critical  buckling  node  for  elastic  general 
instability  failure  of  the  cylinder  with  the  heavy 
frames  replaced  by  typical  frames  (i.e.,  a  cylinder 
equal  in  length,  but  with  no  heavy  frames), 
(see  "List  of  Symbols"  for  other  variables) 
In  computing  PC,  the  radius  to  the  combined  centroid  of  the 
heavy  frame  and  its  effective  length  of  shell  plating,  an 
ELH  had  to  be  determined.   This  was  also  provided  in 
Blumenberg ' s  report : 

ELH=  SL»F1»(AF+SL»T)  +BH 

AFK-f  SL*T 

It  can  be  seen   to  be  a  form  of  simple  ratioing  of  the 

18 
original  EL  formula  developed  in  the  S~P  analysis 

Another  feature  of  HVYFBM  is  to  send  a  value  of  zero  back  to 

the  main  program  if  the  heavy  frame  spacing  computed  by 

GINST  is  equal  to  (or  greater  than)  the  hull  length.   This, 

of  course,  would  mean  that  heavy  frames  are  not  required,  and 

that  the  hull  end  closures  are  used  in  their  place. 
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CHAPTER  NILE 

WEIGHT- DISPLACEMENT. RATIO 

The  simple  weight/displacement  ratio  was  selected  for 
the  optimization  oriterion.   Not  only  was  it  the  simplest 
to  use,  but  it  also  seemed  to  be  the  most  general,  all- 
encompassing  determinant  of  an  optimum  design.   Due  to  the 
subprogram  system  utilized  in  the  optimization,  another 
form  of  criterion  could  easily  be  substituted  if  the  need 
arose. 

Essentially  (see  figure  1?),  WTD3P  computes  the  weight 
of  a  hull  section,  the  length  of  which  is  equal  to  one 
heavy  frame  spacing,  splitting  the  heavy  frame  on  each  end 
in  half.   This  weight  is  divided  by  the  same  hull's 
displacement  in  sea  water,  taking  into  account  internal  or 
external  frames. 
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CHAPTER  TEN 
MAIN  PROGRAM 


Essentially,  the  main  program  performs  vei?y   few 
calculations.   Its  primary  task  is  manipulation  of  the 
various  subprograms  and  iteration  control  of  the  entire 
process  so  as  to  arrive  at  an  optimum  design.   See  figure  23 
for  a  basic,  easier-to-follow  main  flow  diagram.   Figure  2^ 
is  the  complete  main  program  flov?  diagram.   In  addition,  it 
computes  "reasonable"  entering  scantlings  for  the  iteration 
process,  based  on  input  data  and  proven  design  practice. 
The  objective  of  this  type  of  main  program  design  was  to 
obtain  a  system  in  which  the  various  modes  of  failure  and 
safety  factors  could  be  changed  or  interchanged  easily  as 
desired.   As  all  failure  modes  are  separated  into 
subprograms,  it  is  possible  even  to  substitute  an  entirely 
different  analysis  for  any  individual  failure  mode.   Safety 
factors,  which  so  often  are  subject  to  modification  due  to 
new  test  data,  are  also  easily  replaced  or  changed. 

The  input  data  required  (see  Appendix  D  for  format) 
are:   overall  hull  length,  bulkhead  spacing  (ray  or  may  not 
be  given;  if  not,  it  computes  an  optimum  B3),  internal  or 
external  framing  desired,  Young's  modulus-,  metal  yield  stress, 
metal  density,  and  Poisson's  ratio.   In  addition,  as  presently 
set  up,  the  program  may  be  Iterated  for  various  input  values 
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of  depth  and/or  hull  diameter. 

The  depth  input  (in  feet)  is  converted  to  a  collapse 
pressure  PC  (in  psi)  requirement  by  using  the  equation  of  the 
mean  line  drawn  in  figure  18;  this  represents  an  average 
depth  vs.  psi  curve  for  the  various  oceans  included.   The 
figure  (minus  the  mean  line)  was  taken  from  the  handbook  of 
Ooean  Engineering  Tables,  published  by  the  U.  3.  Naval 
Oceanographic  Office,  and  compiled  by  E.  L.  Bialek. 

The  frame  constant  CC  is  computed  from  the  equation  of 
the  curve  shown  in  figure  19*   CC  is  used  to  determine  frame 
proportions  to  be  used  for  various  depths.   The  curve  is  an 
average  of  proportions  of  many  ring  stiffeners  used  in 
present  generation  submersible s.   The  method  of  obtaining 
frame  proportions  is  printer"!  in  the  program  output  (see 
Appendix  D).   A  more  detailed  explanation  may  be  found  in 
Appendix  A. 

The  program  next  computes  its  ''datum  point",  ox- 
starting  set  of  hull  scantlings.   Shell  thickness,  T,  is 
computed  from  the  simple  hoop  stress  formula: 

T=  I'C^pn 

2.0*SIGY 

Frame  spacing,  or,  more  accurately  in  this  case, 

unsupported  length  of  shell  between  frames  (SL)  is  computed 

from  the  parameter  0,  where:     A   */_ , ,   _,,*  '-— - 

9=  ^3(1-?)  ^ 


Figure   1 8 


68 


o 

O 

O 

o 

9 

o 

N" 

(0 

Q- 

UJ 

O 

Ll 
O 

h 

o 
o 

U. 

Id 
C3 

Z 

< 

O 

Ld 
ex 
3> 
<n 

Id 
c^: 
Q_ 


(Id)  HldBQ 


69 


O 

CO 

i 


16 

\ 

1 

F 

1 

RAME  C 

HMENSK 

3M  CUR\ 

/E 

i4 

\ 

i7 

\ 

IO 

\ 

\ 

8 

\..t 
V 

• 

^ 

--— _ 

6 

t 

- 

IOOO 


Zooo 


3O0O 


4ooo 


5  coo 


GOOO 


DEPTH  (FT) 

;    are  19 


4/70 


70 

Reference  6  lists  the  usual  range  of  0  for  present  day 
submersible s  as  1.0-2.5*   r£he  main  program  (see  below)  is 
set  up  so  that  it  starts  at  a  maxmium  frame  spacing  and 
then  decreases  it  in  further  iterations.   Thus,  a  high 
starting  value  of  G  would  be  desirable.  After  several 
calculations  at  various  depths  and  geometries,  a  value  of  0 
=5*0  was  used  for  the  datum  point.   This  gives  a  wide  range 
of  SL  values.   The  iterative  process  of  the  program  produces 
values  of  9  that  eventually  go  low  enough  to  bracket  the 
above  range*   Other  methods  of  obtaining  a  start:'!  rig  SL  were 
investigated,  such  as  the  combined  solution  of  hoop  stress 
and  the  Windenberg-Trilling  equation,  but  all  gave  too  wide  a 
range  (usually  too  large  a  frame  spacing)  as  depth  increased. 
Thus  8  ,  which  does  not  vary  greatly,  was  selected. 

A  starting  frame  size  is  selected  based  on  the  ratio  of- 
frame  cross-sectional  area  to  the  cross-sectional  area  of 
one  unsupported  length  of  shell  plating.   Reference  6  lists 
the  normal  range  of  this  ratio  to  be  0.2-0.8.   During 
rather  extensive  investigation,  however,  it  was  found  that  a 
ratio  of  greater  than  0.5  gave  frames  that  were  grossly 
overdesigned.   Thus,  the  starting  ratio  was  set  at  0.5 • 

AF=0.5*SL*T 
The  frame  moment  of  inertia  is  then  computed  in  accordance 
with  relations  developed  in  Appendix  A.   Originally,  it 
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was  intended  to  select  a  starting  frame  size  by  a  more 
accurate  approach.   However,  the  only  formula  which  could 
be  solved  in  anything  approaching  a  cloyed  form  was  far  too 
conservative  (i.e.,  Tokugawa*s  formula,  reference  25).   The 
advantage  of  a  computerized  approach,  i.e.,  investigation 
of  a  wide  range  of  variables  at  great  speed  and  low  cost, 
was  the  deciding  factor  in  using  the  more  random  iterative 
method  described  above. 

At  this  point  the  program  goes  into  a  double  loop, 
iterating  on  3L  (outer  loop)  and  AF  (inner  loop).   The 
scantlings  are  tested  by  RNLDS.   If  they  do  not  fall  in  the 
lobar  buckling  mode  at  design  collapse  depth,  the  program 
goes  on  to  THKNS.   If  the  scantlings  fall  RNLDS,  SL  is 
decreased  by  multiplying  it  by  0.9,  and  the  RNLDS  test  is 
rerun.   The  decision  to  run  RNLDS  at  this  point  was  made  in 
order  to  start  as  near  the  ''shoulder"  (or  at  least,  not  to 
the  right  of  it)  of  figure  12  as  possible,  since  this  is 
generally  acknowledged  to  be  an  area  of  optimum  design. 
Also,  the  decision  to  vary  SL,  rather  than  T,  to  achieve  a 
set  of  scantlings  that  would  not  fail  the  RNLDS  test,  was 
made  for  two  reasons : 

1)  From  figure  20,  it  can  be  seen  that  either 
T  or  SL  could  be  successfully  used  to  change 
RNLDS  failure  pressure,  with  T  bein 
slightly  more  effective. 


Figure   20 
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2)  The  next  subprogram  in  line  is  THKNS.   It  is 
connected  to  RNLDS  via  another  iteration  loop. 
To  have  both  programs  converging  on  different 
values  of  T  would  cause  endless  loops  in  almost 
every  case.   Thus,  RNLDS  was  iterated  on   SL. 
The  decision  not  to  vary  AF  to  change  RNLDS  failure  pressure 
is  justified  by  figure  21.   The  Nott  buckling  equation  is 
almost  identical  to  RNLDS,  and  thus  can  be  considered 
the  same  for  these  qualitative  investigations.   It  may  be 
noted  that  for  some  geometries,  AF  has  no  effect  at  all  on 
PCR;  thus,  AF  was  not  used  in  the  RNLDS  loop.   This  same 
reasoning  (see  figure  22)  was  used  in  deciding  against  using 
an  AF  variation  of  any  kind  for  convergence  within  the 
THKNS  subprogram. 

As  stated  above,  once  the  scantling  set  (revised  or 
unrevised)  gets  successfully  past  RNLDS,  it  goes  on  to 
THKNS.   There,  the  hull  thickness,  T,  is  adjusted  so  as  to 
have  the  shell  fail  by  9.x i symmetric  yield  exactly  (within  1%) 
at  design  collapse  pressure.   At  that  point,  the  scantlings 
with  revised  T  are  again  looped  back  through  RNLDS  to  insure 
against  lobar  buckling,  and  SL  adjusted  again  as  necessary. 
If  SL  has  to  be  re-adjusted,  the  scantlings  again  go  through 
THKNS.   The  process  continues  until  either: 

1)  a  set  of  the  same  scantlings  pass  without 
change  through  both  RNLDS  and  THKNS,  or 


EFFECT  OF  FRAME  SIZE  ON  ASYMMETRIC 
CKUNG   PRESSURE  BY  NOTT  .  .  . 


7^ 


O.Z 


tte^h 


OA  0.6  0.8 

FJ2AME  AREA/SWELL  AREA 

(NOTE:  TAKEN.  FROM  REFERENCE    15) 
Figure   21 


I.O 


4-/70 


EFFECT  OF  FRAME  SIZE  OiM  AXISYNSMETRIC  75 
(YIELD)  FAILURE  PRESSURE  BY  LUNCIIIK.  .  . 


300  O 


a 
^^ 

ui 

O- 


< 


»ooo 


MAT'U-   MX-  80 

HULL   DIAMETER.  .  16  FT 


F5=  8' 


t.G 


50  io. O 

FRAME  AfcEAOu*) 


15. o 


Figure   22 


J^GlMltY 
4/70 


76 

2)  SL  becomes  so  small  that  there  is  less  than 
four  inches  clearance-  between  adjacent  frame 
flanges. 
In  the  latter  case,  the  loop  is  skipped  without  print-out, 
and  the  next  iteration  is  started,  much  the  same  as  what 
occurs  if  THKNS  cannot  internally  converge  on   a  shell 
thickness  (see  chapter  six). 

Once  the  scantling  set  gets  past  THKNS,  T  is  again 
adjusted  upward,  if  necessary,  by  STRTHK.   A  loop  to  re-test 
back  through  RNLDS  and  THKNS  is  not  used  here,  and  would 
serve  no  useful  purpose,  since  T's  from  THKNS  and  STRTHK 
could  rarely  ever  be  made  to  converge  (i.e.,  the  largest  T 
from  both  programs. is  used). 

Prior  to  entering  GINST,  the  integer  L  is  set  equal  to 
0.   This,  used  as  one  of  the  entering  arguments  for  GINST, 
is  a  control  variable.   3>0  insures  normal  (iterative  BS 
design)  operation  of  GINST.   L=l  is  usely  solely  to  obtain 
a  single-pass  value  of  ii  for  HVYFRM  (see  chapter  eight ), 

The  scantling  set  then  enters  and  is  tested  by  GINST. 
BS  is  adjusted  as  necessary  to  insure  PCG  falls  between 
1.05  PC  and  2PC,  the  only  exception  being  if  BS  becomes 
greater  than  total  hull  length  (see  chapter  seven)*   The 
reason  that  BS,  rather  than  T,  SL,  or  AF,  is  used  here  to 
converge  on  a  satisfactory  PCG,  is  because  in  virtually  any 
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cage,  regardless  of  depth  or  geometrical  proportions,  33 
has  a  much  smaller  effect  on  hull  weight  for  the  same  PCG 
reduction  or  increase.   Since  hull  Weight  already  had.  been 
optimized  for  shell  yield  failure  (the  desired  situation), 
any  further  T,  3L,  or  AF  changes  would  most  probably  take  the 
design  seriously  off  the  optimum. 

Once  past  GIiIST,  the  scantlings  pass  through  KVTFRft 
to  obtain  a  suitable  heavy  frame  design  and  WTDSP  to 
calculate  that  particular  scantling  set's  relative  efficiency 
ratio. 

At  this  point,  values  of  the  scantlings  a.re  printed  out. 
Generally  speaking,  the  print-out  of  each  line  of  scantlings 
can  be  said  to  be  optimized  on  thickness  of  the  shell, 
although  certainly  other  factors  (3L,  BS,  etc.)  change  as 
necessary  to  keep  the  design  on  a  yield  failure  basis. 

dice  out  of  the  first  loop  cycle,  AF  is  decreased  by 
multiplying  by  0.8.   This  is  done  for  ten  cycles,  so  that 
the  original  AF  is  reduced  to  0.108*AF  in  the  last  cycle. 
SL  is  reduced  in  the  same  manner  in  the  outer  loop,  and  as 
described  above,  also  may  be  reduced  within  each  iteration 
as  -necessary  to  pass  RNLDS.   At  the  end  of  both  loops  (100 
major  iterations)  the  program  terminates  for  that  set  of  input 
parameters,  and  prints  out  the  optimum  (i*e.,  smallest) 
weight/di splacement  ratio . 


('BEAD 

1  DATA 


COMPUTE  INITIAL 
SCANTLING  SET 


11 

<J=1,10 

IE 

9 
S=1.10 


78 


PCR=RNLDS 
(ASY.BKLG. ) 


CALL 


. THKNS 
(AXI.YLD. ) 
1 


-9 

\        STH'PHK 


G    I  L 


GINST 
(GEN. INST. ) 


HVYFRM 


WD=tfTDSP 


*: 


3U 


WRITE     I 

ut;  utJ 


REDUCE 
FRAME 


SL 
IEDUCED 


SIMPLIFIED 

MAIN 
PROGRAM 


REDUCl 
SL 


$ 


are   23 


Alio 


0<i/  18/7C 
CHART    TlILe    -    FFCCECURES 


*UTCFICS  CHSRI  SET  - 


E.S.pcginley:    suaiiulopt  *aiN         taGE      02 


I 


01 


/    KEAij    FP.OM    OEV    / 
/  5  / 

/       VIA    FC"AI  / 

•  133  / 

/    IMC    Tl-c    LIST    / 


I 
I 


•  USI    ■    KIUM 

•  tU.    f.    £,    SI 

•  «!-.,",      .-.VJ 

•  ••»«••      i 

1 

• 

1 

03 

/     Rt»U    FROM 
/                  5 
/       lit    FORMAT 
/                   '00 
/    INTO     THE    LIST 

CEV   / 
/ 
/ 
/ 

1 

LIST    =    N, 

(OPTH  [(ItN)iNN 

1,M 


/    READ    FROM    OE V    / 
/  5  / 

/       Via    FCRMAT  / 

/  20C  / 

/    INTC    THE   LIST    / 


1 

I       V.OTE       06 

44*44444 

*    *    *• 

•               LIST    =    M, 

• 

*           1  OP  II  PM  I,  MM 

=           » 

«                   1 ,  M  > 

* 

***    4    s    4    r    4 

1 

*     *     * 

1 

07 

/ 

/ 

/    WRITE     TC    CEV       / 

/             t 

/ 

/       III    FORMAT 

1 

/                    l« 

1 

/    WRITE    TC    DEV       / 


/     via   fcr»at 

/ 

/ 

/ 

/ 

FRCP    THE    LIST 

/ 

1 
1 

1      NCE 

09 

* 

4*44*44 

* 

*    * 

• 

LIST     =     R  HO .     S I GY 

f    * 

* 

6,    GNU 

* 

» 

4      *      4      *      *      *      * 
1 
1 
1 

*    * 
10 

♦ 


•      11 

BS    .EQ.     0.0       * 


I  12 

* * 

\         6SF    ~    BS/W.O         I 


1 

13 

/    WRITE    TO 
/                    * 

/     via  format 

/                 850 
/    FRCM    THE    LIST 

CEV 
/ 

/ 
/ 

t 

1 

1 

1       *CT£       14 

•  LIST    =     ESF  4 

**4***44444 
I 
>l 

20  I       F.CTE       15 

444444444*4 

4       CCNTI'I'JE        4 

44444444444 

I 

I        kCTE       16 

•  4444444444 

•  BEGIN   00    LOOP  4 

•  125    II    -   li    k        4 


04.  ?7 >  | 

I 


17 


I  DP1H    =    OPTHHI  I  I  | 

I  I 

I  CC    =  I 

I  1C.0«EXP(  (  320.0    -  I 

I  CPThl/1000.01    ♦  I 

I  6.5  I 

« 4 

I 
I 

I  18 

/    WRITE     TC    OEV       / 
/  6  / 

/      via    f CPPAT         / 
/  7  50  / 

/     FRCM    THE    I  1ST    / 


1 
1 

MOTE 

19 

4       « 

» 

4      4      4 

4      4      4 

«      • 

* 

LIST    = 

CC 

• 

1 

NOTE 

zc 

BEGIN   OC    LCCP 
225    JJ    =    1,    H 


I 


21 


I  BS   =    BSO                 I 

I  I 

I  OM    =    DMK  JJI            I 

I  I 

|  RP    -    DP/2.0             I 

I  I 

|  OKF    =    DM/12. 0          I 


I 


Figure   2h  a 


•         22 

•         • 

*                    • 

t                              • 

TRUE 

^   .ne.   t.o 

4-4 

1 

•              * 

| 

*      * 

1 

• 

.FALSE 

3    . 

' 

01    . 

i 

16 

t 

i 

23 

/    WRITE    TO    OEV       / 
1  b  I 

I     via  FOPMai       / 

/       250       / 
/  FROM  THE  LIST  / 


4  LIST  •  0 


NOTE   24 
4  4  4  4  4 

THt  DMF   * 


04/18/ 7C 
CHAPT     I  I  iLt     -     PFCCECUKES 


AUTCFLCh    CHAP!     St  I 


00 

E.S.1CG1NLEY:    SuriHULOPT    main  PAGE      03 


/    «S|TE 
/                    6 

/     vu  ru.- 

/                  35C 

/  fcch   in:  i 

10 

■■I 

CEV 

/ 
/ 

/ 
/ 

1 

1 

1 
1 

r.j 

It 

02 

•    IIS1 

■    ->FT 

'• 

3«F 

• 

02.24- 
17 

1 

—  >l 
1 

0} 

/  mizc  ia  o;v     / 


/     vi»  f  dukai 

/  *50 


l 


0- 


I  PC    =  I 

I  O.CC5E75»UPTM»'-2/  | 

I  1J40C.O    >  | 

I  0.445»DPTh  | 

I  | 

I  T    =  | 

I  PC»D>/(2.P«SI(,Y1  I 

I  I 

I  SL    =  | 

1  J.O'SCKUKMI)  | 

I  /( n.o»<  i.o  -  I 

I  CNU«»2I l«<0.25l  | 


OS 


I  AF    =    0.5«SL»F  I 

I  t 

I  B    »  I 

I  SURIIAF/12. 7«CC I  I  | 


I 


1                   CCDP    ■ 

1 

1          <2.89*CC»«4    • 

1 

1          1  1.34>CCs»3    » 

| 

1            1?.47»CC"2     ♦ 

1 

1                 *.36«CC1 

1 

1    /(U.O'CC    «    20. 

•1      1 

1       Fl     =    CCCP>8««4 

1                     1    «    I 

1 

1 

1 

07 

1                Tl    -    I 

1 

I          WDCPT     =    130.0 

1 

1        KCIt 

OH 

•    •*•*»*>« 

*    * 

•          BEGII.    DC    LCOP 

» 

•          1 1    J    ■     1,     1C 

* 

•    **«**«*« 

*  * 

04.21 >l 

1       NOTE 

09 

•    ***«     ***» 

«     * 

•          BEGIN    DO    LCOP 

* 

•             t    K     •     1,     10 

• 

•     *•     ««...( 

*    • 

04.18  — ->l 

1 

10 

I              AFM(K)     =    Af 

1 

1              S1MJ!     ■     SL 

1 

I                 FU    -     CC>6 

1 

1 

1 

>«< 

I       ^IE      11 
•    »*«»••• 

CCM1MJE  • 

I 


•  LI.     IFk 
4.01 


I 


FNLCSIT.SL  I 
I 
I 


14 


THK^ 

(I.SL.FI.PC, 

CIS! 


I     ALTERNATE     PETUP.N        I 

* * 

I        IS  4.15     I 


I  Si    =    0.")«SL  | 


SL     .GT.     IFW    «     • 
»  4.01  • 


I 
I 
I 
I 
I  20 

» 9 

I  SL     =     SL/0.9  | 


Figure   ?J\  '■ 


04/18/70 
CHART     Hill     -     PROCICUStS 


AUTCFLCk    CHA-M     !ET 


81 


E.S.KCGINLEY:    Subkuiopt    »A|N  net      04 


0J.I5-  — >| 


Sibihk 
ll.Sl.f  I, PCI 


1 

11 

1 

C  IF-S1 

1. 

11,51, Fl.ec, 

10 

F>CO.NI>F,,BS,L  , 

11 

KULKTH 

11 

H 

13 

fuvfrm 

H 

1. 

(T.SL.FI, 

H 

10 

HLLMH.FCPC& 

h 

11 

BS.AFF .8F.AF, 

H 

1 

FCH 

H 

1 

1 

1                         C6 

kOO.JI     *  I 

WIOSCI I, SI ,FI,BS,  I 

MI.Ki,."    I  j 

WCC   ■    kC(K,J)  I 


«  * 


»  X    .CI .    I 


I 


I       W00PT    •    kEll.l)       I 


I 


C4.07* — >* 
15  • 


WDIK.JI     .LT. 
»         WDOPT 


I       HOCPI    -    kCIK.J)        | 


-04.09* — >| 
14  I 


I  FS    =     SL     ♦     B  I 

I  I 

I  BSF    =    BS/12.0  I 


I 


13 


/    plP.1  IE     ID    DEV       / 
/  6  / 

/      VIA    FORMAT         / 
/  55C  / 

/    FRCM    THE    LIST    / 


I       NOTE       14 


LI  SI    ■    I,    FS,    e, 
BSF,     Bll,     KUB 


03.12«-->l 


K    .IT.     10 


IFAL  SF 

I 


u 


I             SL    =    SI  Ml  1  I  I 

• • 

I 


13  I  I 

I    ■    1 

AF     =     0.8'»FM(K| 

e  « 

SQ'OISF  /<2.7*CCI  ) 

fi  =  cccp«a»«4 


5  •       is 


•  END    CF     OC 

•        LOOP? 


si  --  o.e«SLi(ji 
t  .   n 

»F    =    0.5*SL*1 

B    • 
SOKIUF/I  2.7»CC>  I 


I       FI    =    CC0r«3«»4  | 


11  •         21 


•         ENC    OF    DO 
•  LOOP? 


l»ES         .       3    . 


•  •       TRUE 
WDOPT    ,NE.  • 

•  1C0.0         • 

•  * 

*         * 

I FAL  SF 


/    WRITE    TO    OEV      / 


/       VIA    FORMAT  / 

/  6  49  / 


IB  |  ?4 

/    WRI TE    TC    OEV      / 
/  6  / 

/       VIA    FORMAT  / 

/  650  / 

/    FROM    THE    I  1ST    / 


I       BOTE       25 
LIST    =    WOCPT  • 

I 


•  END    DF     DO 

»     .   LOOP? 


125  •  27 


ENO    CF     DO         «-» 

LOOP?  «  | 


F3    ure   24   c 


CALL    Ex  IT 
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CHAPTER  ELEVEN 

RESULTS  AliD  CONCLUSIONS 

Various  parametric  studies  were  performed  with  the 
program,  once  it  was  checked  against  some  contemporary 
submersible  hulls  to  see  if  it  indeed  was  "in  the  ball 
park".   The  most  obvious  study  was  to  see  how  the  W/D 
ratio  varied  with  depth  (figure  25).   This  was  done  for 
three  common  steels  in  use  today,  all  with  approximately 
"ideally  plastic"  stress-strain  characteristics.   The  points 
obtained  plotted  into  smooth  curves  on  the  semi-log  plot 
used.   It  is  fairly  obvious  that  if  a  W/D  greater  than  0.5 
is  considered  unsatisfactory,  the  following  limits  would 

appiy : 

STE2L  •  LBi^^CJ^Jf/p_^L  0 

HTS  3500  FT. 

HY-80  6200  FT. 

HY-100  7^00  FT. 

A  second  study  was  conducted  to  observe  how  the  W/D 
ratio  varied  with  hull  diameter,  all  other  factors 
(including  depth)  being  constant.   It  was  found  that  it 
varied  very  little  if  at  all,  as  can  be  seen  from  the 
following  set  of  data: 
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DEPTH  .THOUSANDS  OF  FEET 
Figure   25 
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Given:-  HY  80  steel,  collapse  depth=5000  ft.: 
HULL^piAJ-lBTlgR  (FT.)  W/D .PATIO 

20  0.417 

30  0.418 

50  0.41? 

80  0.409 

It  could  tentatively  be  said  then,  that  for  a  given  depth, 
hmll  efficiency  (W/D)  will  he  approximately  the  same 
regardless  of  how  large  the  hull  diameter  may  be. 

Figure  26  shows  the  result  of  plotting  the  variance  of 
W/D  with  metal  yield  stress.   This  also  plotted  into  a 
smooth  curve  in  this  case.   It  is  probable,  however,  that 
for  higher  yield  strength  (and  therefore,  strain-hardening) 
metals  requiring  a  slightly  different  analysis,  the  curve 
would  have  a  sharp  break. 

The  third  plot  attempted  (figure  27)  at  first  appeared 
to  be  a  hopeless  scattering  of  data,  but  after  some  analysis 
revealed  rather  interesting  results.   The  points  plotted 
were  taken  from  a  single  optimization  (one  diameter,  one 
depth)  such  as  in  the  example  print-out  in  Appendix  D.   A 
card  computing  thinness  factor  was  inserted  in  the  program 
and  printed  out  with  the  regular  data. 

The  ,!SL  loop"  is  actually  a  series  of  M/'D   points 
computed  for  scantling  sets  with  the  same  fj   ■  i   >-  c i ng , 


OPTIMUM  WT/A  RATIO  VS.   STEEL  STRENGTH 
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CONSTANTS  • 
COLLAPSE     DEPTH-      IOOO  FT 
HULL     DIAMETER-  3o  FT 

HULL     LENGTH:        200  FT 
FRAMES  :    INTERNAL 
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YIELD    STRESS  THOUSANDS  OF  PSl 


Figure   26 
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Figure  2? 
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buy  varying  frame  size  as  an  input.  In  each  case,  it  is 
strongly  evident  that  a  saddle  point,  or  optimum  W/D  for 
that  frame  spacing,  was  reached.  This  means  that  for  at 
least  this  particular  set  of  inputs,  the  program's  frame 
size  (AF)  iteration  range  was  large  enough  to  bracket 
optimum  values.  One  of  these  frame  spacing  miniraums,  then, 
was  the  optimum  W/D  ratio. 

The  "W/D  envelope"  encloses  all  W/D  values  computed  in 
the  program.   There  seems  to  be  good  indication  that  the 
"optimum  W/D"  indicated  is  a  true  optimum,  since  the  lower 
portion  of  the  envelope  rises  on  either  side  of  it. 

Another  indication  given  by  the  plot  is  that  the 
program  gives  a  much  wider  variation  in  W/D  ratios  with 
larger  frame  spacings,  for  varying  frame  sizes.   As  frame 
spacing  becomes  smaller,  the  W/D.  ratios  produced  become  more 
"convergent". 

Perhaps  the  most  obvious  conclusion  from  figure  26  is 
that  A  alone  is  certainly  not  an  accurate  indicator  of 
optimum  W/D  ratio,  although  it  could  be  utilized  (after 
extensive  data  gathering)  to  indicate  the  general  area  in 
which  to  design. 

It  is  interesting  to  note  that  average  computation  • 
time  for  each  optimum  computed  was  less  than  one-half 
minute.   Also,  it  was  interesting  to  note  that  the 
optimum  solutions  contained  as  scantlings,  generally 
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speaking,  smaller  frames,  larger  heavy  frame  spacing,  and 
slightly  thicker  shells  than  submersibles  of  present  design. 
Evidently,  this  results  in  better  V//D  ratios. 

One  of  the  major  advantages  of  this  program  and  its 
general  method  of  computation  and  print-out  (see  Appendix  C) 
is  that  it  gives  a  great  variety  of  alternatives  from  which 
to  choose.   For  instance,  perhaps  the  optimum  V//D  ratio 
for  a  certain  hull  configuration  and  depth  contained 
scantlings  which  gave  a  very  small  frame  spacing. 
Analytically,  this  may  give  a  superior  W/D  ratio.   But 
practically  speaking,  the  cost  of  fabricating  and  installing 
a  great  number  of  frames  may  be  out  of  the  question.   Thus 
(particularly  if  the  submersible  is  not  critically  weight- 
limited),  the  table  of  printouts  may  be  "browsed"  for 
more  attractive  scantling  sets:   ones  with  acceptlble  W/D 
ratios,  but  with  inherently  lower  construction  costs.   It 
should  be  repeated  here  that  each  of  these  printed  lines 
are  not  mere  random  choices.   Each  line,  prior  to  print-out, 
has  already  been  through  one  of  the  main  hurdles  of  the 
optimization  program.   The  design  of  shell  thickness  and 
frame  spacing  combination  to  give  failure  by  axi symmetric 
yield.,  and  thus  most  efficient  use  of  the  material's 
inherent  strength,  is  completed  prior  to  printing  each  line 
of  the  answer  table.   Thus,  the  pro -'ram  may  be  used  not  only 
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as  a  structural'  optimization,  but  also  is  an  indi sponsible 
reference  for  any  economic  analysis  of  a  particular  design. 
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CHAPTER  TWELVE 
RECOMMENDATIONS 


This  Program. 


!•      Fianp:esr_of  Variatlorie      One   of   the  most   common 

questions  after  a  particular  run  completion  was 
"I  wonder  if  this  is  rj§allv_  the  optimum,  or 
didn't  I  go  far  enough  in  frame  (or  frame  spacing) 
variation?"   It  is  "believed  that  plots  such  as 
figure  27  for  each  run  could  probably  give  a 
definitive  ansx-jer  in  most  cases e   Mot  only  can 
one  tell  if  each  set  of  frames  ran  through  an 
optimum,  but  it  should  also  usually  be  possible 
to  tell  if  the  "W/D  envelope'*  passed  through  its 
optimum.   This  may  be  a  rather  painful  way  to 
assure  one's  self  that  his  run  covered  all  the 
territory  that  was  necessary,  but  at  this  writing, 
it  appears  to  be  the  surest  way.   As  an  alternative, 
some  rather  extensive  studies  could  be  conducted 
using  various  depth  and  other  inputs  to  test  range 
validities.   Once  determined,  the  program's  range 
of  AF  and  SL  variation,  governed  by  the  indices  K 
and  J,  respectively,  could  be  controlled 
appropriately.   From  a  rather  cursory  inspection  by 
the  author,  it  appeared  that  in  runs  conducted 


91 


in  data-gathering  for  this  paper,  the  ranges  of 
J=K~10  were  adequate  in  all  "but  a  few  questionable 
cases. 
2»   Effects  of  Fabrication  Procedures*   It  is  rather 
obvious  that  without  any  doubt,  the  greatest 
shortcoming  in  the  optimization  is  in  its 
allowances  for  such  relative  "unknowns"  as  residual 
welding  stresses,  low-cycle  fatigue,  shell  or 
frame  out-of -roundness,  etc.   The  gross 
approximations  made  by  REDPR  and  STRTHK  certainly 
fall  far  short  of  the  accuracy  of  the  various 
failure  mode  analyses.   There  is  no  doubt  that  more 
basic  research  must  be  accomplished  in  this  region 
before  a  completely  dependable  optimization 
program  could  be  achieved.   As  it  stands  now,  the 
safety  factors  built  into  the  program  could  fall 
into  two  categories: 

a)  Completely  safe  design,  in  which  all 
scantlings  are  overdesigncd  to  the  extent 
that  design  "optimization"  is  almost  useless. 

b)  more  realistic  (i.e.,  lower)  safety  factors 
based  on  scanty  experimental  data,  which 

is  not  universally  applicable,  and  thus 
might  be  considered  unsafe. 
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It  is  hoped,  that  the  safety  devices  employed  in 
this  program  adhere  to  a  "middle  o?   the  road." 
policy.   REDPR,  it  is  believed,  results  in 
somewhat  of  an  overdesign  for  the  two  between- 
frame  failure  modes,  but  in  underdesign  for  the 
general  instability  mode.   This  underdesign  is 
hopefully  picked  up  in  STRTHK,  which,  as  mentioned 
in  Chapter  6,  insures  "a  reasonable  stress  level 
in  the  frame  flange  prior  to  collapse  for  those 
frames  with  an  initial  out-of -roundness"1  . 

The  only  way,  it  appears,  to  resolve  these 
questions,  is  in  extensive  testing  of  models  with 
deliberate,  measurable  defects  of  a31  types.   The 
most  needed  data  of  this  type  is  in  the  failure 

mode  most  effected  by  defects:   general  instability. 

p 

It  is  believed  by  some  authors   that  the  out  of 

roundness  analyses  developed  thus  far  are  overly 
pessimisltlc  when  applied  to  full  sized  submersibles. 

The  most  obvious  extension  of  this  program  is  into 
strain-hardening  materials.   Eecause  of  the 
program's  general,  characteristics  (i.e.,  a  main 
flow  control  nrogram  manipulating  subprograms 
which  actually  do  pressure  calculations),  it  can 
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be  easily  transformed  with  -very  few  changes  other 
than  subprogram  additions. 

The  need  for  such  a  transformation  is  obvious, 
when  one  observes  hull  steels  projected  for  future 
(and  in  many  cases,  present)  usage.   Any  steel 
with  a  yield  strength  greater  than  100,000  psi 
can  be  considered  to  be  a  strain-hardening  material. 
Details  of  the  method  of  transforming  the  present 
program  into  one  suitable  for  use  with  high-strength 
steels  may  be  found  in  Appendix  E. 
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APPENDIX  A 
FRAME  DI MEM S I ON 3 
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APPENDIX  A 

TYPICAL  FKAHE 3: 

In  order  to  arrive  at  reasonable  frame  cross-sections, 
averages  of  a  great  variety  of  frame  scantlings  used  in  a 
number  of  submersibles  at  various  depths  were  obtained  (see 
figure  19 t  chapter  10).   The  common  denominator  for  all  of 
these  frames  is  the  web  thickness,  B.   This  method  was  taken 
from  Adamchak  in  reference  1 ,  who  did  the  same  thing  for 
surface  ship  frames. 

The  ratios  used  in  figure  19  vjere  taken  from  the 
following,  derived  from  the  averages  computed  (see  figure 
28): 

Flange  width=FW=CC*B,  where  CC=y  of  figure  19. 
Flang e  dep th=FD=F W=  COB 
Flange  thickness=1.7*B 
In  addition  to  the  above,  the  remaining  terms  in 
figure  28  are  defined  as  follows:. 

FC:    Distance  of  frame  centroid  from  plate  neutral 

axis 
DN:    Distance  of  combined  centroid  from  plate 
neutral  axis 
The  various  formulas  for  frame  and  combined  frame-plate 
its  of  inertia  '  id  for  center  of  gyration  were  tal  sn 


T-  SECT!  ON 
DIMENSION     IDENTIFICATION 
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from  a  general  development  of  these  in  reference  19.   They 
are  listed  in  easily  programmed  form,  and  according  to 
reference  I9i  differences  from  values  of  actual  standard 
production  frames  due  to  fillets,  etc.,  are  always  less  than 
two  per  cent. 

The  combination  of  these  two  methods  enables  the 
computer  to  easily  arrive  at  any  variety  of  T-stiffener 
characteristics,  each  of  which  is  proportioned  according  to 
present  submersible  design  practice. 

After  substituting  the  expressions  for  FD,  FW,  TF  and  B 
into  the  equations  for  T-frames  in  reference  19»  the 
following  expressions  result  for  typical  (i.e.,  not  heavy) 
frames,  and  are  used  in  the  program: 
Frame  area,  AF; 


AF=FW*TF+FD*B  =  2.7*CC#B**2 
Frame f mmome^t_pfmi  inertia^ ^FI^ 
FI    =    CCOP#B*#^,    where: 

CCOP=  (2 .  89*CC#*4+11  ♦  3^*CC**3+13.47*CC**2 
+^.36*CC)/(12oO«-CC+20.^) 
Effective  frame-plate  moment  of  inertia  (using; 
effective  length  of  shell  plating; ,  EL,  developed 
in  reference  18). -El: 
EI=  0 .  225*CC*B***H*  ( CC+1 .  7  )  **2#  ( ( 1 .  07*CC+0 .  562 )/ 

(  CCH  1 .  7  )+3.  0-  ( 1 .  6"3+(T-l .  7*B  j/(  (  CC+1 .  7  )  ■:. B  )  ) « *2/ 
(1.0+(2.7*-CC<  •  2)/(E]  -.  r))) 
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FC: 

FOB-;;-  (  ( CC+1 .  7  3/1 .  227-1 .  0/1 . 1?6  )+'J?/2.  0 
DN: 

DN=FC*2 . 7*CC*B**2/ ( 2 . 7*CC*B**2+EL*T ) 
HEAVY  FRAMES: 

Heavy  frames  were  averaged  in  the  same  manner  as 
above,  although  less  data  exists.   The  relationships 
obtained  did  not  depend  upon  depth,  and  were  determined  to 
be  as  follows: 

FDH=17,0*BH 
FWH=13.0*BH 
TEH=2.0*BH 

cco?=i 330.0 

Heavy  frame  area, _AV'd: _ 

AFH=i|-3.0*BH**2 
^BXI  ,£rg:-?.e.  fflPIl1^.^.  °f .— ipertia.^  FIII.; 

Filial  330.  0-s:-}3-::-«-'l- 
FCHj.. 

FGH«1^«  25#BH+T/2 . 0 

DNH=FCH*AFH/(AFH+ELH*T) ,  where 

ELH-eff ective  length  of  shell  plate  next  to  the 
heavy  frame  (see  chapter  8) 
EIH  was  not  need  eel  in  the  program. 
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APPENDIX  B 
LUNCHIK'S  PLASTIC  KINGS  ANALYSIS 
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APPENDIX  B 
LP/JGHIKJ.S  MASTIC  HINGE  _AI'ALYSIS 
(Refei-  to  reference  lh   for  this  discussion) 
Lunchik^s  final  Pc/Py  equation  is  developed  through  the 
use  of  his  parameters  B  and  K.   There  are  not  solvable  except 
by  assumptions  which  Lunchik  makes.   Two  of  his  assumptions 
(i.e.,  kpyukx  and   /c^K/Kp)   are  easily  worked  out  from 
substitutions  in  identities  from  his  paper.   However,  one  is 
less  clear  (i.e.,    ^/d^^     ^)»  particularly  because  B^ 

is  misprinted  as  fy   in  reference  1^. 

Ik 

From  Lunchik1 s  definitions   in  his  analysis  of  a  one 

unit  square  element: 

(l)M0=kjj!p  =  circumferential  edge  moment 
(2)  |SL~ Kdph E  circumferential  compressive  membrane 
force  • 

(4)oyrA  =  |<.pE  circumferential  membrane  stress 
Assuming  (i.e.,  approximating)  the  circumferential 
bending  stress  to  be  elastic: 

^>b0-  T>;/  bk3^  2   bh*   (DH1    ^ 
From  (2)  and  (k) : 

Using  (5)  and  ( 6) : 

3S}    hz   Np   h-K^pk   n  Kj$ 
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To  further  re-arrange  these  assumptions  to  obtain 
expressions  for  Bx  and  B0  to  use  in  program  ELNCK: 
(7)  fy-i^'/d^  ;  (8)  \C*/\i,-<f«*/d^ 

b^fk,  nt  "TfK ^  ^      (9) 

(7)i  (Q)»  and  (9)  are  then  put  into  Lunchik's  equations 
(12),  (13),  and  (15)  for  9, ,  ©a  and  %   ,  and  these  in  tux- 
are  substituted  in  his  (25)  and  (26),  which  are  used  to 
solve  for  the  plastic  reserve  stength  ratio  Pc/Py 
(PCTR  in  program  ELNCK). 
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APPENDIX  C 


PROGRAM  LISTING 
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RING-STIFFENED  CYLINDRICAL  SHELL  OPTIMIZATION 


.. .  INPUT  PARAMETERS... 
MATERIAL  DENS  I TY  =  0 , 2 B 5  LBS/CU. IN.   YIELD  ST  RES S=l OOOOO,  PSI 
YOUNGS  MODULUS^  0.3E  06  PSI     POISSONS  RATIO=0.3o 


*FRAME  DIMENSIONS  MAY  BE  OBTAINED  BY  USING  THE  FRAME 
CONSTANT  BIWE8  THICKNESS)  AS  FOLLOWS: 
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•••LIGHT(TYPlCAL)  FRAMES... 
FRAME  DEPTH= FLANGE  WIDTH=   6.6  X  B 
FLANGE  THICKNESS^  1.7  X  B 

...HEAVY  FRAMES... 
FRAME  DFPTH=  17  X  B 
FLANGE  WIDTH=  13  X  B 
FLANGF  THICKNESS^  2  X  B 


DEPTH=    5000.0    FT.       DIAMETER=    30,0    FT, 


FRAMES: INTERNAL 


.. .DESIGN 

OPTIONS... 

TYPICAL 

TYP  ICAL 

HEAVY 

HEAVY 

WEIGHT/ 

SHELL 

F  R  A  M  E 

FR^ME 

FRAME 

FRAME 

DISPLACEMENT 

HICKNFSS 

spacing 

CONSTANT-: 

«   SPACING 

CONSTANT* 

RATIO 

(  IN.  ) 

(  IN.  ) 

(  IN.) 

(FT.) 

(IN.  ) 

5.0350 

108.09 

3.438 

200,00 

0.000 

0.565 

4.  5805 

107.73 

3.075 

200.00 

0.000 

0.501 

4.  5324 

107.40 

2.7  50 

200.00 

0.000 

0.474 

4.2713 

107.11 

2.460 

200.00 

0.000 

0.435 

4.0395 

106.85 

2c  200 

198.16 

2.981 

0.411 

3.8703 

106.62 

1.968 

73.77 

2.6  64 

0.399 

3,  6163 

106.41 

1.760 

47.17 

2.463 

0.  378 

3.6163 

106.2  2 

1,  574 

38.32 

2.354 

0.374 

3.6163 

106.06 

1.408 

?0*64 

2*21  5 

0.393 

3- 5768 

105.91 

1.260 

20.64 

2.084 

0.378 

4.6269 

86.80 

3.075 

200,00 

0.000 

0.531 

4.3938 

86.47 

2.750 

200.00 

0.000 

0.484 

12? 

4„ 196? 

86,19 

2.460 

200.00 

0.000 

0 .  4 '.  6 

4. 1139 

85.92 

2.200 

194.23 

3.082 

0.433 

3.7526 

85.69 

1.968 

79.  98 

2.677 

0.401 

3.7247 

85.48 

1.760 

65.73 

2.651 

0.390 

3.5618 

85.29 

1.5  74 

44,41 

2.399 

0.374 

3.5618 

85.  13 

1,408 

30.22 

2o  285 

0.  378 

3.5228 

84.99 

1.260 

23.14 

2.154 

0.  376 

3.5228 

84.  85 

1,  127 

16.07 

2.014 

0.  385 

4.2901 

69.73 

2.750 

200,00 

0.000 

0.505 

4.0739 

69,44 

2,460 

200,00 

0.000 

0.460 

3. 8^08 

69.18 

2,200 

200.00 

0.000 

0.422 

3°  72  66 

68,94 

1.9  68 

190.39 

2.956 

0.401 

3. 5705 

68.74 

1.  760 

70.10 

2.652 

0.388 

3.4143 

6  8  o  5  5 

1,574 

58.67 

2,514 

0.365 

3. 4143 

68.39 

1.408 

35.99 

2.354 

0,368 

3.4143 

68.24 

1,2  60 

24.51 

2.  199 

0.372 

3*4143 

69.10 

1,127 

18.83 

2.074 

0.374 

3.5850 

67.98 

1.009 

13.17 

1.999 

0,402 

4.0094 

56.04 

2.460 

200.00 

0.000 

0.482 

3.8607 

55.7  8 

2.  2  00 

200.00 

0.000 

0.444 

3.6871 

55,55 

1.968 

199.97 

3.068 

0.416 

3»5315 

5  5,34 

1.760 

80a06 

2.660 

0.  396 

3.  2  61.1 

55.  16 

1.574 

66.27 

2.596 

0.  362 

3.3121 

54.  99 

1.4C8 

43.36 

2.385 

0,  362 

3,4397 

54.84 

1.260 

34.22 

2c320 

0.  369 

3.4536 

54.71 

1.127 

25.10 

2.153 

0,370 

3.4677 

54,59 

1.008 

16.01 

2.020 

0.383 

3.5216 

54.4  8 

0,901 

11.47 

1.971 

0.404 

3c 6047 

45,06 

2.200 

94.08 

3.014 

0.474 

3.4426 

44,  93 

1.968 

94.08 

2.845 

0.430 

3. 3388 

44.62 

1.760 

86.65 

2*  684 

0.  398 

3.0832 

44.44 

1.574 

68.13 

2.590 

0.  362 

3.2879 

44.2  7 

1.408 

64.44 

2.537 

0.364 

3.3768 

44,12 

1.2  60 

38.70 

2.335 

0,369 

3. 3905 

43,99 

1.127 

31,37 

2.222 

0.365 

3.4822 

43.87 

1.008 

20.40 

2.097 

0.380 

3.4963 

43.77 

0.901 

16.75 

1.971 

0,379 

3.5105 

43.67 

0.8  06 

13.12 

1.921 

0.388 

3.22  94 

36.26 

1.968 

67  .  5o 

2.849 

0.449 

2.9829 

36,05 

1.760 

67.  50 

2.708 

0.398 

2.9176 

3  5,87 

1.5  74 

67.50 

2.603 

0.367 

3.  1114 

35.70 

1.4C9 

67.50 

2.549 

0.363 

3. 19  55 

3  5.55 

1  .  2  60 

43,  80 

2.359 

0.  362 

3. 2819 

3  5.42 

1.127 

34.95 

2.286 

0,  362 

3.3706 

3  5,30 

1,0  03 

26,12 

2.  130 

0.366 

3.3943 

35-19 

0.901 

17.33 

2.002 

0.  374 

3.4757 

35.10 

0.  8C6 

14.40 

1.928 

0,382 

3. 48^9 

35-01 

0,721 

11.48 

1,902 

0,  392 

2.8816 

29.19 

1.760 

60.14 

2.755 

0.420 

2.8272 

29,01 

1.574 

60.14 

?o  661 

0,385 

2.9036 
3.  0154 
3.0969 
3. 1806 
3. 2666 
3.4439 
3.3788 
3.4313 
2.6059 

2. 6770 
2.8109 

2.8869 
3.0312 
3.1131 
3.  1973 
3.2837 
3.3725 
3.4249 
2.34  98 
2. 5830 
2.7121 
2.8477 
2.9901 
3. 1396 
3.2245 
3.3117 
3.4012 
3.4150 
2.3059 
2.  5  72  5 
2.7011 
2.8362 
2. 9780 
3.1269 
3.2114 
3. 2982 
3.3873 
3.4011 


28.  84 
28.  69 
2  0.56 
28.  A 4 
28.33 
28,24 
28.  15 
28.08 
23.52 
23.35 
23.  21 
23.07 
22.95 
22.85 
22.7  5 
22.67 
22.59 
22.52 
18.97 
18.82 
18.68 
18.56 
18.46 
18,36 
18.28 
18.20 
18.13 
18.07 
15.31 
15.  17 
15.05 
14.95 
14.85 
14.77 
14.69 
14.62 
14.56 
14.51 


1,4  08 


1. 

2  60 

1. 

127 

1. 

0  08 

0. 

901 

0. 

806 

0. 

721 

o. 

645 

1. 

574 

1, 

400 

1. 

260 

1. 

127 

1. 

008 

0. 

901 

0. 

8  06 

0. 

721 

0. 

645 

0. 

5  77 

1. 

408 

1. 

260 

1. 

127 

1. 

008 

0. 

9C! 

Oo 

8  06 

0. 

721 

0. 

645 

0. 

5  77 

0. 

516 

1. 

2  60 

1. 

127 

1. 

008 

0, 

901 

0. 

806 

0. 

721 

0, 

645 

0, 

5  77 

0. 

516 

0. 

461 

60.14 
55.36 
38.  70 
29.22 
22.13 
17,43 
12.73 
12.73 
55.86 
55.  86 
55.86 
40.47 
34.74 
25.22 
17.63 
13.85 
11.97 
11.97 
37.26 
37.26 
37.26 
35.71 
28,02 
21.90 
15.81 
14.29 
12,78 
11.27 
31.68 
31.68 
31.68 
31.68 
24.25 
18.10 
14.43 
13.21 
12.00 
12.00 


2.600 
2.441 
2.282 
2.  171 
2-071 
1.949 
1-890 
1.848 
2.683 
2.626 
2.488 
2.301 
2.238 
2.086 
1.969 
1.894 
1.865 
1.832 
2.531 
2.42  9 
2.327 
2.246 
2.130 
2.044 
1.914 
1.855 
1.830 
1.821 
2.425 
2.  357 
2.283 
2.  195 
2.059 
1.954 
1.874 
1.833 
1.817 
1.789 


128 
0.370 
0.360 
0.  358 
0.360 
0.365 
0.3  78 
0.382 
0.  378 
0.399 
0,375 
0.360 
0,356 
0.355 
0,357 
0,367 
0.376 
0,384 
0.381 
0.  389 
0.376 
0.  360 
0.  355 
0.356 
0,  363 
0.371 
0.373 
0,379 
0.383 
0.  387 
0.382 
0.366 
0.355 
0.357 
0.  366 
0.372 
0.  374 
0.380 
0.374 


OPTIMUM    WEIGHT/DISPLACEMENT    RATI  0:0, 355 
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APPENDIX  E 

METHOD  OF  PROGRAM  CONVERSION 
TO  ACCOMMODATE?  STRAIN-HARDENING  MATERIALS 
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APPENDIX  E 
OUTLINE  OF  METHOD  FOR  CONVERTIMJ  TO  AN 

OPTIMIZATION  PROGRAM.  FOR  STRAIN-HARDENING.  J^IATERIALS 

Strain-hardening  metals  differ  from  elastic  perfectly 
plastic  metals,  in  that  above  the  yield  point,  the  stress  is 
not  a  constant  value  as  strain  increases,  but  increases 
(usually  at  progressively  slower  rates)  as  strain  increases. 
This  means  that  true  plastic  flow  is  never  achieved;  the 
metal  continues  to  retain  a  modulus  of  some  value,  albeit 
smaller  than  the  original  constant  value.   This  results  in  a 
type  of  combination  elastic-plastic  buckling  failure  in 
stiffened  cylindrical  shells,  termed  inelasti_c  failure. 
It  is  particularly  important  to  analyze  subersibles 

constructed  of  strain-hardening  materials  with  strain- 

11  13  20 

hardening  analyses   *  •"  '    .   Normally,  a  submarine  hull  is 

designed  to  have  some  plastic  yielding  (i.e.,  initially 
beginning  in  the  hull  adjacent  to  the  frame  flanges) 
somewhere  between  operating  and  collapse  depth;  indeed,  there 
must  be  yielding  prior  to  collapse  depth  in  order  for  the 
hull  to  crush  there.   With  ideally  plastic  materials,  once 
the  yielding  point  is  reached,  generally  speaking,  buckling 
failure  is  ruled  out  (or  it  would  have  occurred  earlier,  due 
to  hull  g e o m e t ry ) .   W i  t h  s t ra in- ha r d e n i ng  ma t e r ia 1 s , 
hex;  per,  c  tc  -  tl   hull  b  gins  bo  pL  i  bically  deform, 
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buckling  failure  is  not  ruled  out.   In  fa^ct,  there  is  a 
strong  possibility  that  the  buckling  failure  (which  will 
occur)  will  happen  at  a  lower  pressure  than  the  yield  failure 
calculated  by  an  ideally  plastic  analysis,  due  to  i/ne 
reduced  metal  modulus*   Thus,  to  use  an  .ideally  plastic 
collapse  pressure  analysis  on  a  strain  hardening  material 
might  give  dangerously  overoptimistic  failure  predictions, 
particularly  if  the  metal  has  a  very  high  yield  point. 

For  convenience  in  analyses,  strain  hardening  materials1 
stress-strain  curves  are  characterized  not  only  by  E,  but 
also  by  E-fc  (ETAN,  tangent  Modulus)  and  Es  (ESEC,  secant 
modulus).   See  figure  32. 

Because  with  this  type  of  stress-strain  curve,  the 
modiili  are  always  dependant  upon  the  stress  state  when 
above  yield  stress,  it  is  necessary  for  all  straiii-hardening 
collapse  pressures  to  be  computed  using  an  iterative  process „ 

It  is  generally  the  approach  to  the  solution  of  this  process, 

20 

and  an  example  pressure  analysis  by  Reynolds   that  will 

comprise  the  rest  of  Appendix  E. 

Essentially,  the  critical  collapse  pressure  is  obtained 
when  the  buckling  equation  (depending  upon  buckling  node, 
references  11,  13,  or  20)  is  solved  simultaneously  with  the 
pre-buckling  equation  (stress  intensity  as  a  function  of 
prnrs-r.  ,-: )  ,  where  : 


tfc 


Figure   32 
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Stress  intensity^  6-^\6l*6l -6X  6S      (see  figure  33) 
Either  V~G  or  S-P  stress  theories  could  be  used  to  calculate 
the  stress  intensity,  S.   V~G  is  less  accurate,  but  3-P 
might  give  convergence  problems,  depending  upon  hull 
geometry  and  depth  at  which  stress  is  calculated. 

One  difficulty,  that  of  finding  a  way  of  describing  a 
strain-hardening  stress-strain  curve  with  a  minimum  of  input 
data,  is  solved  in  reference  23.   in  this  method,  the  entire 
curve  may  be  approximated  with  extreme  accuracy  by  using  only 
four  inputs  (see  figure  3*0  :   E,  dy,  0a  and  ob  .   By 
manipulating  some  of  the  Romberg-Osgood  equations,  it  is 
possible  to  obtain  E^  and  Es,  given  any  value  of  stress,  via 
the  following  expressions: 

loSlo(<l/tfb) 
E/Et=  l+0.^2857n(d/c5'Jn*', 

E/Eg=   1+0.^2857(6/^)^ 

This  will  be  assumed  to  be  the  content  of  a  subprogram 
called  Subroutine  ROHOS  (see  figure  35). 

Both  the  pre -buckling  equilibrium  equation  and  the 
buckling  equation  will  be  approximated  in  the  region  of 
interest  with  straight  lines.   This,  and  some  of  the  following 
methods  of  determining  the ' intersection  of  the  two  equations, 
were  patterned  generally  after  similar  methods  used  in 
reference  22. 
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METHOD  OF  OBTAINING  ROMBERG -OSGOOD 
INPUT  PARAMETERS  FOR  COMPUTING  Es  AND  Et 


Figure    "}h 
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COMMON 

/G/   SIC 

rA,SIGB| 

EN=1.0+0, 

-3853/A1 

.OGIOCSIGA/SIGB) 

1. 


ENl=SN-l-.0 
I 


ETAN=E/  ( 1 .  O-JO .  42857*EN*  ( SIG/SIGA )  *#EN1 ) 


ESEC=E/ ( 1 . 0+0 . ^2857# ( SIG/SIGA ) #*EN1 ) 


FLOW  DIAGRAM: 
SUBROUTINE  ROMOS 


Figure  35 
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Subroutine  LINE.   This  subroutine  computes  EMM,  the 
line  slope,  and  BEE,  the  line  stress  intercept  of  the  pre- 
buckling  equation,  using  the  S-P  analysis  (see  figure  36). 
The  values  given  by  LINE  will  change  only_  with  scantling 
changes.   This  subprogram  would  constitute  the  only  deviation 
of  the  main  flow  program  from  the  perfectly  plastic  case. 
It  would  be  placed  in  the  main  program  directly  after  point  3 
(i.e.,  prior  to  the  RNLDS  call).   It  would  also  be  used 
within  THKNS  whenever  T  changes. 

The  following  discussion  will  involve  the  solution  of 
the  asymmetric  inelastic  buckling  mode  as  developed  by 
Reynolds  in  reference  20.   The  sane  general  iterative  process 
would  be  utilized  in  the  solution  of  the  other  two  modes  of 
hull  failure. 

Function  PCRP,   See  figure  37 •   This  small  subprogram 
merely  computes  asymmetric  inelastic  buckling  pressure, 
using  as  inputs  ETA?!  and  ESEC  (computed  by  ROKOS),  T,  SL, 
and  PCREl  (computed  as  PCRE  by  the  elastic  portion  of  PCRP's 
calling  program,  RNLDS).   The  pressure  is  computed  using 
Reynolds'  equations  outlined  in  reference  20  (for  a 
simplified  presentation,  however,  see  reference  12,  which 
gives  the  axisymmetric  mode  also). 

Subroutine  RNPT.   See  figure  38.   This  subprogram  obtains 
the  intersection  of  the  pre-buckling  equilibrium  equation 
with  the  line  determined  by  the  two  input  pressures  (PRSl 
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COMMON/D/E , GNU , DM , Z , CC , SIGY , CCOP , RHO/H/BEE , EMM 
.  '/U/THETA  ,  ALFA ,  BETA __ 


I 


P=0.9*PC 


± 


CALL 
FRAME 
V(FI,T,SL) 


I 


10 

=1,2 


SIGU=-P*RM/T 


CALL 
PULOS 
(SL.T.FI.P.Fl 

^,F3,F4,A) 


SIGJ  iP=SIGUv;-  ( 1 . 0-A*F2  ) 


SIGi-lX=0.5*SIGU 


S ( I ) -SQRT ( SIGMP#*2+SIGMX**2-SIGMP#SIGMX ) 


I 


P(I)=P 


d^D 


YES 


->- 


P=PC 


no 


EM=(S(2)-S(1))/(P(2)-P(1)) 


I 


BEE~-S(l)-EHi^:-P(l)l 


FLOW  DIAGRAM:  SUBROUTINE  LINE 


Figure  "}6 


COMMON  /D/E , GNU , DM , Z , CC , SIGY , 
CCOP.BHO  


3 


GNUP=0 . 5- ( 0  c  5-GNU ) *ESEC/E 


I 


PCRP= ( PCRE1* ( 1 . 0-GNU**2 ) *ETAN/ ( 1 . 0-GNUP**2 ) ) 
( 1 . 0+( 0 . 922*SQRT ( RM*T ) /SL ) * ( ESEC/ETAN-1 . 0 ) ) 
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FLOW  DIAGRAM: 
FUNCTION  PCRP 


Figure  3? 
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COHMOiJ/D/E  ,GNU ,  DM ,  Z  ,  CC  ,  SIGY  ,  CCOP ,  RHO ] 
/G/SXGA , SIGB/ii/BEE  .EMM 


± 


CALL 


ROIIOS 
(3IG2.ETAN, 
ESEC) 


i 


PRS2-PCRP ( PCRE1 , T , ESEC , ETAN , SL ) 


1 


ENEMSIG2-SIG1  )/(PRS2-PRSl ) 

IGEE^SIG  1  -ENNftPRSl 


P= ( CEE-BEE ) / ( EMM-ENN ) 


FLOW  DIAGRAM: 
SUBROUTINE  RNPT 


Figure   38 


47-70 


140 


and  PRS2)  and  stress  intensities  (31G1  and  SIG2)  in  terras 
of  the  intersection  coordinates,  S  and  P.   The  values  EMM 
and  BSE,  which  determine  the  pre -buckling  equation,  are 
piped  into  RNPT  via  a  COMMON  statement,  along  with  other 
normal  program  data,  inputs,  including  SIGA  and  SIGB.   Input 
arguments  include  SlGl ,  S1G2  and  PRS1 ,  which  along  with  the 
inelastic  buckling  pressure  at  point  2,  PRS2  (computed  by_ 
RNPT  from  input  PCRE1),  describe  end  points  of  the  buckling 
.equation  line  approximation. 

Function  RNLDS.   This  subprogram  is  the  same  as  RNLDS 
used  in  the  ideally  plastic  case  (i.e.,  computation  of  elastic 
lobar  buckling  pressure,  PCREl ) ,  with  the  addition  of  a 
programmed  iteration  (figure  39 »  as  used  in  reference  22)  for 
the  inelastic  portion.   The  first  two  stresses  used  in  the 
iteration  are  SIGY  and  0.9---3IGY.   When  the  difference  of 
computed  inelastic  failure  pressure  and  assumed  stress  at 
which  failure  will  occur  in  the  next  iteration  becomes  less 
than  one  half  of  one  per  cent  of  the  present  predicted 
failure  pressure,  convergence  is  assumed.   The  iteration 
sequence  may  be  followed  by  using  the  diagram  of  the  inelastic 
portion  of  RNLDS  (figure  40)  with  the  plot  (figure  39). 

As  noted  before,  identical  procedures  would  be  followed 
in  the  axisymmetric  and  general  instability  cases.  The  main 
progra  i  would  be  unalte   ."   it'  the  exception,    ■  ■  ■  as  noted 


ITERATIVE  METHOD  USED  TO  CONVERGE 

ON  INELASTIC  COLLAPSE  PRESSURE  .  . 
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(CONTINUING  ON  FROM 
PRESENT  RNLDS  PROGRAM) 


VlZ 


COMMON   /G/   SIGA.SIGB   /H/  BEE,    EHH  | 


KcG'MLtY 

Alio 


I 


Sl=0.9-»-SIGYl 

t 

CALL 

ROIIOS 

v(Sl,ETAN,ESECj/ 

I 


Pl-PCRP  ( PCRE1 ,  T ,  ESEC  ,  JSTAN  ,  SL ) 


I 


S2=SIGY 


I 


CAJ  L 
RNPT 
v(PCREl,T,SL,Pl,Sl,S2f> 

P2.P3.S3) 

—   * 
CALL 

RNPT 

\|;pcre1,t,sl>p2is2,/ 
S3.f^,P5.s4) 


c 


P3-P^  >  0  i  005.:;-P3^iP ^JPCRP1~P4 


v  YES 


CALL 
RNPT 

V(PCRE1,T,SLIP^,S3, 

S^,P6,P7,S5) 


^P5-p6>o.oo5^p5^°— . >-[pcrpi~?6 


♦  YES 
CALL 
RNPT 
V(PCREltTlSLfP6,S^, 
S5,P8,P9,S6) 


=P8 


/P7-P8  >  0 .  0Q5-:--p"8>VNO    >-|PCRPl^P 
|  YES 


CALL 
RNPT 
V(PCRE1,T,SL,P8,S5 
S6.P10,P11,S7) 


e 


j-  3  Mire  40 


y^s 


P9-P10>0. 005*P9?VM->4PCRP1=P10 


(•CONTINUED  FROM 
IAST  PAGE . . . ) 

i 


1^3 


WON- CONVERGENCE 


RNLDS=REDPR(PCRE1 tPCRPl ) 


(Figure  ^0  Continued) 
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before,  of  insertion  of  subprograi-i  LINE.   It  is  obvious 
that  an  optimization  using,  the  inelastic  analysis  would  take 
much  longer  (perhaps  by  five  or  six  times)  than  the  ideally 
plastic  case. 
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